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SECTION  1 

GENERAL  INFORMATION 


1.1  -OB;JECTIVE 

The  purpose  of  this  document  is  to  provide  information  about 
the  Special  Operations  Radio  Frequency  Management  System  CSORFMS) 
model  enhancements  produced  by  Systems  Exploration,  Inc.  The 
enhancements  include  the  incorporation  of  a  Sounder  Update  CSU) 
model  and  an  improved  version  of  the  Field  Strength  CFS)  model. 
The  SU  model  is  an  adaptation  and  translation  to  FORTRAN  of  the 
Army  PROPHET  Evaluation  System  (APES)  BASIC  model.  The  FS 
enhancements  involve  modification  of  the  field  strength  modeling 
algorithms,  consequential  to  efforts  to  improve  the  accuracy  of 
this  model . 

SORFMS  will  be  run  on  a  Grid  Compass  computer  or  any  of  a 
family  of  IBM  compatibles.  The  SORFMS  upgrades  were  developed, 
integrated  and  tested  on  an  HP  9050  computer.  Since  the  source 
code  is  FORTRAN  77,  it  can  be  transported  to  the  hardware 
environment  with  little  or  no  modification.  However,  since  there 
are  no  facilities  for  comprehensive  testing  in  this  environment, 
it  is  the  HP  9050  version  that  is  described  herein.  The  HP  9050 
software  is  referred  to  as  High  Frequency  Tactical  Decision  Aids 
(HFTDA)  or  simply  TDA. 

The  methodology  of  using  the  HP  9050  as  a  developmental  tool 
for  software  destined  for  the  Grid  Compass  and  IBM  compatibles  is 
modeled  after  a  similar  method  used  for  developing  Advanced 
PROPHET  (AP).  Similarities  among  SORFMS,  AP  and  TDA  have  been 
exploited  wherever  applicable. 

This  document  is  addressed  to  analysts  and  maintenance 
programmers  who  must  understand  the  rationale  and  details  of 
these  enhancements,  exercise  the  functions,  and  integrate  these 
new  features  with  additional  changes  to  the  programs  produced  at 
NOSC.  For  continuity  among  the  relevant  documents  (listed  in 
1.5),  Section  1  herein  provides  an  overview  of  SORFMS, 

1 . 2  BACKGROUND 

The  SORFMS  is  a  small,  lightweight,  stand  alone  and  highly 
transportable  real  time  propagation  assessment  and  forecasting 
system  which  defines  natural  propagation  constraints  on  HF 
transmissions  and  outputs  this  data  in  an  easily  interpreted 
format.  The  system  includes  an  HF  transmitter,  HF  receiver, 
spectrum  monitor,  frequency  management  terminal,  and  a  portable 
computer.  The  system  software,  described  herein,  is  installed  on 
the  portable  computer 

The  SORFMS  software  is  a  collection  of  computer  simulation 
models  developed  to  support  tactical  use  of  the  HF  band  (2  -  32 
MHz).  Primarily,  SORFMS  can  determine  the  existence  of  an  HF 
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skywave  channel  between  two  sites  anywhere  in  the  world;  it  can 
also  determine  the  potential  for  a  hostile  force  to  intercept  the 
transmission,  radio  locate  the  transmitter,  or  jam  the  reception 
site.  For  short-range  circuits,  a  groundwave  model  is  included; 
this  along  with  ^Irhe-  skywave  model,  can  be  exploited  to  provide 
(for  example)  covert  communications  between  closely  operating 
units. 


In  SORFtiS,  most  adjustable  parameters  can  be  preset  and 
will  be  accessed  automatically  as  needed.  There  are  a  few 
parameters  that  must  be  continually  adjusted,  and  these  are  the 
sol ar /geophysi cal  constants  that  drive  the  computer  models.  The 
Sounder  Update  CSU)  model  is  used  to  update  a  pseudo-sunspot 
number  and  is  included  in  this  category.  For  normal  operation,  * 

these  parameters  are  input  in  near  real  time.  However,  for 
planning  purposes  or  to  analyze  past  performuince  on  a  circuit, 
estimated  or  historical  data  can  be  used. 

Reference  1,  paragraph  1.2.1  provides  a  detailed  summary  of 
the  physics  of  HF  propagation.  This  particular  reference 
describes  Advanced  PROPHET;  however,  the  physical  description  is 
directly  applicable  to  SORFMS. 

1 . 3  SUMMARY 

The  enhancements  to  the  SORFMS  described  in  Section  2 
i nc 1 ude: 

a.  Addition  of  a  SU  model.  The  model  was  derived  from  APES  and 
"ranslated  from  BASIC  to  FORTRAN,  the  »SU»  command  was  added  to 
the  AP  command  menu  so  that  it  is  accessible  to  the  operator. 

b.  Modification  of  the  FS  model  to  provide  more  accurate  data 
outputs.  Specifically,  the  FS  model  now  provides  outputs  which 
can  be  more  closely  correlated  with  the  CCIR  data  base. 

1.4  ENVIRONMENT 
1.4.1  Hardware 

SOF^FMS  is  targeted  to  operate  on  the  GRID  Compass  computer, 

IBM  PC  and  compatibles.  The  SORFMS  modules  described  herein  were  ^ 

developed  on  the  HP  9050  computer. 

A  version  of  Advanced  PROPHET  (called  TDA)  runs  on  the  HP 
8050  system  at  NOSC.  This  system  was  used  as  a  test  bed  for 
development  and  test  of  the  SORFMS  SU  and  FS  modules.  The 
adaptation  of  the  8050  program  to  the  GRID  Compass/IBM  PC  and 
compatibles  environment  entails  the  tranferrance  of  the  FORTRAN 
source  code  to  a  machine  which  can  compile  into  GRID  Compass 
machine  code. 
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1.4.2 


Software 


SORFMS  is  configured  to  run  under  the  liS-DQS  operating 
system,  Version  2.0,  or  greater.  The  TDA  development  system  runs 
under  the  UNIX  operating  sy^em. 

1 . 5  REFERENCES 

1.  NOSC  TD  848,  Operational  Users  Manual  for  Advanced  PROPHET  on 
MS-DOS-BASED  Microcomputer  Systems,  by  R.  W.  LaBahn,  August  1985. 

2.  NOSC  TD  1065,  Operational  User’s  Manual  for  Tactical  Decision 
Aids  for  HF  Communication,  by  D.B.  Sailors,  July  1987. 

3.  Integrated  Logistic  Support  Plan  (ILSP)  for  Special  Operation 
Radio  Frequency  Management  System  (SORFMS)  (Draft),  prepared  by 
U.S.  Army  Communications  Electronics  Command,  COMM/ADP  Center, 
Fort  Monmouth,  NJ,  September  1986. 

4.  NOSC  TD  782,  Tactical  Decision  Aids  for  HF  Communication,  by 
D.B.  Sailors,  15  November  1984. 

5  NOSC  Contractor  Report,  Program  Package  Document  for  Sounder 
Update  and  Field  Strength  Software  Modifications  for  Special 
Operations  Radio  Frequency  Management  System  (SORFMS),  prepared  for 
D.  B.  Sailors,  by  J.  R.  Gnessin,  D.  J.  Brandon,  and  D.  L.  Lucas, 
Systems  Exploration,  Inc.,  March  30,  1988.  (NOSC  TD  1848,  Volume 
2) 
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SECTION  2 

PROGRAM  DESCRIPTIONS 


2. 1  OVERVIEW 


The  Sounder  Update  CSU)  and  Field  Strength  CFS)  modules  of 
SORFMS,  described  in  this  section,  are  separate  modules  in  the 
sense  that  each  is  independently  invoked  by  the  operator. 
However,  most  of  the  modules  of  SORFMS  are  interrelated  in  the 
sense  that  they  use  common  routines,  and  in  some  cases,  the 
outputs  of  one  module  may  affect  the  operation  of  another, 
depending  upon  the  sequence  selected.  The  case  in  point  is  SU, 
where  the  sunspot  number  produced  by  the  SU  module  is  one  of  the 
parameters  used  by  the  FS  module. 

The  SU  and  FS  modules  invoke  many  standard  subroutines  and 
library  routines  previously  existing  in  TDA  which  were  not 
modified  and  are  thus  not  described  in  detail  here,  other  than 
references  when  they  are  used.  Subroutines  DAMBOLT  and  MUF85 
were  modified  slightly  to  accommodate  the  new  LTLFLD  algorithm. 
However,  changes  did  not  affect  the  algorithms,  functionality,  or 
compatibility  with  the  remainder  of  the  modules  with  which  they 
interface.  Therefore,  these  subroutines  are  not  described  here. 
However,  for  configuration  management  purposes,  their  source  code 
listings  are  included  with  the  listings  in  the  SORFMS  Program 
Package  Document  (reference  6,  listed  in  paragraph  1.5). 


2.2  FIELD  STRENGTH  CFS) 

2.2.1  Scope 

2.2. 1.1  Identification 
PROGRAM  FSTREN 

2.2. 1.2  Engineer 

D.  L.  Lucas 

2. 2. 1 . 3  Purpose 


This  program,  and  its  associated  subroutines,  provides  a 
method  for  the  rapid  assessment  of  HF  skywave  signal  quality  as  a 
function  of  frequency  within  the  propagation  bandwidth  for  links 
between  specified  points  on  the  earth's  surface. 
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2.2. 1.4  Background 

The  computer  module  development  to  replace  "DAMBOLT"  long 
path  model  for  paths  shorter  than  7000  KM  is  a  model  relying  on 
propagation  via  the  regular  E-layer  (refererTce  1)  and  the  F2- 
layer  (reference  2).  The  routine  finds  a  lower  order  mode  for 
the  E-layer,  r2-layer,  and  a  mixed  E-layer  and  F2-layer 
mode.  Only  one  E-layer  hop  is  considered  in  the  mixed  mode.  The 
remaining  hops  are  F2-layer,  either  below  or  above  the  MUF  for 
the  path. 

The  MUF  model  used  is  that  of  HFTDA  (MINIMUF-85)  for  the  F2 
layer.  The  MUF  for  the  E-layer  is  calculated  as  in  lONCAP 
(reference  3).  E-layer  MUFs  at  night  cannot  be  lower  than  700 
KHz  (reference  1).  The  ionospheric  absorption  equation  is  the 
one  used  in  lONCAP  with  the  near  specular  reflection  losses 
calculated  for  E-layer  modes  at  low  frequencies.  The  over-the- 
MUF  losses  (reference  4)  are  calculated  using  the  Phillips  method 
for  values  of  the  Freq/MUF  ratio  of  approx.  1.4  to  approx.  1.5 
depending  on  ground  distance,  with  values  greater  being 
considered  in  the  ‘‘scatter"  region.  The  losses  above  these 
values  are  estimated  as  in  lONCAP  which  was  taken  from  some  work 
of  JTAC  (references  5,6).  This  means  the  over-the-MUF  losses  are 
a  function  of  season,  sunspot,  geomagnetic  latitude,  local  time, 
and  path  length  since  the  distribution  of  the  MUFs  used  are  a 
function  of  these  variables. 

The  value  of  the  Absorption  Index  I  at  night  has  a  value 
which  is  a  function  of  solar  activity,  which  was  taken  from  the 
work  of  Wakai  (reference  7).  The  value  perviously  used  in  IQNCAP 
was  from  earlier  data  at  high  solar  activity  (I  .ge.  .1).  This 
meant  the  absorption  was  constant  at  nighttime  hours  for  a  given 
operating  frequency  on  a  given  path  and  solar  activity. 

The  minimum-hop  predicted  for  the  E  layer  is  calculated  for 
an  h'E  of  110  KM  and  the  path  length.  The  radiation  angle 
associated  with  this  mode  must  have  an  angle  greater  than  some 
predetermined  minimum.  If  not  dictated  by  the  user,  the  angle  is 
assumed  to  be  at  the  horizon,  i.e.,  zero  degrees. 

The  minimum  F2-layer  hop  must  be  above  some  pre-set  value 
for  the  take-off  angle  as  for  the  E-layer  and  also  penetrate  the 
E  layer  at  the  F2-layer  angle  calculated  using  h'F2  as  calculated 
in  HFTDA.  The  F2-layer  hops  are  increased  until  one  satisfies  the 
above  restrictions. 

The  mixed  mode  is  calculated  using  the  E-layer  height  (110 
KM)  and  the  h'F2  to  determine  a  take-off  angle.  The  frequency 
must  be  supported  by  the  E-layer  at  one  end  of  the  path  and 
penetrate  at  the  other  end  for  this  mode  to  be  possible.  One  E- 
layer  hop  is  permitted  and  the  path  completed  with  F2-layer  hops 
depending  on  path  length. 
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Only  one  minimum  hop  is  chosen  for  each  mode  of  propagation 
since  the  antennas  are  only  chosen  by  operating  frequency  and 
path  length  making  discrimination  between  a  1  hop  and  2  hop  path 
not  practical  until  more  realistic  antenna  patterns  are  used. 
The  value  of  the  transmission  loss  used  to  predict  the  fields  is 
the  loss  associated  with  the  mode  of  least  loss,  i.e.,  they  are 
assumed  to  be  uncorrelated. 

All  losses  calculated  are  a  function  of  geomagnetic 
latitude,  length  of  path,  time  of  day  and  season  since  the  system 
loss  tables  of  lONCAP  (.'reference  3)  are  used  in  the  module  along 
with  the  distribution  of  the  MUFs  (reference  3).  The  routine 
LTLFLD  is  highly  modular  allowing  changes  in  the  predicted 
parameter  with  some  ease. 

The  ease  in  changing  methods  of  calculating  these  parameters 
also  allowed  for  changes  in  over-the-MUF  (reference  8),  low 
frequency  short-paths  and  radiation  angles  to  improve  the  RMS 
prediction  error. 

2.2. 1.5  Restrictions  and  Limitations 

The  module  LTLFLD  is  built  around  analytical  methods  and 
empirical  data  in  the  range  of  three  to  thirty  MHs .  Any 
calculations  or  predictions  outside  this  range  should  be 
considered  highly  suspect. 

Sunspot  numbers  over  150  should  be  used  with  caution  since 
the  critical  frequencies  of  the  F2-'layer  (FoF2)  are  more  related 
to  where  you  are  within  a  given  cycle  rather  than  the  absolute 
sunspot  number.  In  other  words  depending  on  cycle  the  FoF2  may 
not  increase  from  sunspot  150  to  say  200.  Some  scientists 
believe  a  saturation  effect  takes  place. 

The  predictions  are  hourly-medians  of  the  monthly  medians 
correlated  with  monthly  median  running  average  sunspot  number 
(RASSN)  and  much  care  should  be  taken  when  predicting  for  a  given 
day  or  a  few  days  ahead.  Field  strength  predictions  are  weakly 
associated  with  daily  values  of  solar  activity. 

Solar  flux  data  which  are  used  as  an  indicator  of  ionisation 
in  the  F2-layer  may  or  may  not  effect  in  the  same  manner  the  E- 
and  D-layer  which  controls  most  non-devi at i ve  absorption 
calculated  by  this  method.  Correlation  coefficients  are  not 
available  for  the  D-,  E-,  and  F2-layer  to  warrant  daily 
predictions  of  the  field  strength.  The  accuracy  of  the 
prediction  routine;  therefore,  depends  entirely  upon  its  intended 
use,  i.e.,  frequency  assignment  in  the  long  term,  siting,  antenna 
selection,  day-to-day  frequency  use,  absolute  signal 
determination,  etc. 

Care  should  always  be  taken  to  not  use  the  tool  beyond  its 
intended  use. 


2.2. 1.&  Language 

HP  FORTRAN/77 


2.2. 1.7  Computer  Configuration 

HP  9000  Series  500  HP-UX  (HP  9050) 

2.2. 1.8  Models  Referenced 
None 

2.2.2  Numeric  Method 

The  method  for  predicting  field  strength  at  the  receiving 
location  using  module  LTLFLD  is  the  solution  of  "Norton’s" 
transmission  loss  equation  (reference  9). 

The  system  loss  equation  in  terms  of  resulting  field 
strength  is; 

Ef  =  107.2  +  Pt  +  Qt  +  201 ogf  -  Lb f  -  Li  -  Lm  -  Lg  -  Lh 

where;  Lbf  =  32.45  +  201ogF  +  201ogP' 
f  =  operating  frequency  (MHz) 
pf  =  virtual  slant  range  (KM)  (reference  1) 

Pt  =  effective  radiated  power  in  same  units  as  received 
power  (Watts) 

Li  =  ionospheric  absorption  loss  below  MUF  (dB)  (reference  1) 
Lg  =  ground  reflection  losses  at  intermediate  reflection 
points  (dB)  (taken  as  2  dB  per  ground  reflection) 

Gt  =  antenna  power  gain  relative  to  isotrope  in  free  space  (dB 
Lm  =  over-the-MUF  loss  (dB)  (reference  4) 

Lh  =  excess  system  loss  to  allow  for  auroral,  sporadic-E 
obscuration  and  other  losses  not  explicitly  included  in 
the  predictions  (dB)  (reference  1) 

The  above  losses  are  calculated  separately  and  are  unique 
allowing  any  specific  loss  calculation  to  be  changed  without 
affecting  the  others.  LTLFLD  is  merely  a  solution  to  the  above 
equat i on . 

The  ionospheric  absorption  equation  used  is  that  of  ITSA-1 
(reference  1)  revised  to  account  for  nearly  specular  reflection 
at  low  frequencies  in  the  day-time  hours  as  in  Radar-C  (reference 
4) . 

The  E-layer  critical  frequencies  (FoE)  are  calculated  using 
the  absorption  index  "I"  as  in  ITSA-1. 

The  simplist  calculations  were  used  to  increase  speed  of 
calculation  as  accuracy  was  assumed  not  to  suffer  for  medion 
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values  and  field  strength  calculations  alone.  Fourier  series 
expansions  representing  the  predicted  ionospheric  coefficients 
are  not  included  in  any  of  the  variables. 

2.2.3  Module  Design  Description 

2.2.3. 1  Required  Input 
User  entries; 

1.  Date  of  model  run  Cset  by  the  TDA  command  'DA')  and 
stored  in  variables  MONTH,  JDAY,  YEAR. 

2.  Sun  spot  number  CSSN)  of  model  run  (set  by  the  TDA 
command  'SS'). 

3.  Transmitter  site  parameters  for  the  site  (e.g.,  hono  for 

Honolulu)  including  the  site  latitude  (STALAT CK+X ) ) ,  longitude 
(STALON (K+X ) ) ,  antenna  type  (STAANT CK+X ) ) ,  and  power 

CSTAPWR (K+X ) .  (Values  set  by  the  TDA  command  'AS'). 

4.  Receiver  site  parameters  for  the  site  (e.g.,sdiego  for 

San  Diego)  including  the  site  latituae  (STALAT (K+X )) ,  longitude 
<STALON(K'+X)  ) ,  antenna  type  (STAANT (K+X )) ,  and  power 

(STAPWR (K+X ) ) .  (Values  set  by  the  TDA  command  'AS'). 

5.  Atmospheric  Noise  (ATMOS)  (set  by  the  TDA  command  'AT'). 

6.  Si  anal -to-noi se  ratio  (SIGNSE)  (set  by  the  TDA  command 
'SN'). 

2. 2. 3. 2  Processing 

The  Field  Strength  module  consists  of  a  program  and 
subroutines.  The  program  is  invoked  via  the  executive  routine  in 
TDA  upon  entry  o'f  the  FS  command  by  the  operator.  On  completion 
of  processing,  the  TDA  is  invoked  once  again  and  control  is 
returned  to  the  TDA  executive. 

Figure  1  shows  both  the  hierarchical  structure  and  calling 
sequence  of  the  routines  with  the  Field  Strength  module.  Arrows 
point  to  the  called  routine  from  the  calling  routine.  Each  call 
is  not  indicated,  and  in  general  does  not  follow  the  same 
pattern. 

2. 2. 3. 3  Output 

The  output  of  the  Field  Strength  module  is  a  table  (and/or 
plot)  of  the  field  strength  in  dB  for  the  specified  path  and 
date.  Field  strength  is  provided  as  a  function  of  of  frequency 
and  time.  Frequencies  range  from  2  to  30  MHs  and  time  is  over  a 
24  hour  period. 
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OVERALL  CALLING  STRUCTURE 


ARROWS  POINT  TO  THE  CALLED  ROUTINE  FROM  THE  CALLING  ROUTINE. 
EACH  CALL  INCLUDES  ARGUMENTS.  DATA  COMMUNICATED  THROUGH 
COMMON  IS  NOT  INDICATED.  AND  IN  GENERAL  DOES  NOT  FOLLOW  THE 
SAME  PATTERN. 


Figure  1.  Field  Strength  Module  Hierarchy  and  Calling  Sequence 
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2.2.4  Subprogram  Design  Descriptions 

2.2.4. 1  Subroutine  LTFLD 

2.2.4. 1.1  LTLFLD  inputs 

Call  variables 

FREQ  current  operating  frequency  (MHz) 

CPNT  control  point  array  as  defined  in  MUF85  (deg  fx  KM) 

Common  variables 

PWRWAT  transmitter  power  (watts) 

SF2MUF  F2-layer  MUF  from  MUFB5  (MHz) 

IFHOPS  minimum  F-layer  hops 

HPF2  virtual  height  of  F-layer  from  MUF85  (KM) 

MISC(l)  current  month 

2.2.4. 1.2  LTLFLD  Processing 

a.  Figure  2  illustrates  the  flow  of  LTLFLD.  Following  the 
declarations,  data  is  initialized  and  local  variables  assume 
values  of  corresponding  global  variables.  Subroutines  ELAYER, 
FLAYER,  and  XSLOSS  are  called  in  sequence  to  compute  E  layer  MUF 
(EMUF),  absorption  index  (ABSI),  and  critical  frequency  (F0F2). 

b.  If  the  frequency  is  less  than  1.4*EMUF,  some  E  layer 
reflection  is  assumed  and  subroutine  MODLOS  is  called  to  compute 
the  transmission  losses  associated  with  this  mode.  If  the 
frequency  is  equal  to  or  greater  than  1.4*EMUF,  then  no  E  mode 
reflection  is  computed  for  these  conditions.  After  calling 
MODLOS,  LTLFLD  corrects  for  partial  penetration  of  the  E  layer 
and  totals  the  transmission  losses;  field  strength  is  obtained 
for  this  mode. 

c.  FREQ  is  tested  against  the  frequency  necessary  to 
penetrate  the  E  layer  (E  penetration  frequency  is 
FOEMIN*SECD (2) *1 . 05) .  If  FREQ  is  high  enough  to  get  through, 
calculate  an  F  mode  by  calling  MODLOS  with  F  arguments; 
otherwise,  FHOPS  is  increased  and  a  new  geometry  is  obtained  via 
GEOHOPS  to  try  again.  Up  to  three  tries  are  permitted. 

d.  If  FHOPS  >  1  a  mixed  mode  (lE-NF)  propagation  path  can 
exist.  Assume  total  hops  equal  to  FHOPS,  one  of  which  will  be  an 
E  hop.  Calculate  a  weighted  average  height  and  get  remaining 
geometry  from  GEOHOP.  Test  that  the  ray  path  will  penetrate  the 
E  layer  at  one  end  (F  mode  result),  but  not  the  other  (E  mode). 
If  true,  a  mixed  mode  exists;  calculate  its  characteristics  via 
MODLOS.  I f  no  such  mode  exists,  go  to  last  step. 

e.  Finally,  LTLFLD  sets  the  predicted  field  strength  to  the 
greatest  of  those  determined  thus  far.  The  last  thing  to  do  is 
check  the  calculated  E-layer  MUF  against  the  F-layer  MUF  from 
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c 


LTLFLD 


DECLARATIONS 

AND 

INITIALIZATIONS 


PWRDB  «  LOG(PWRWAT)*  4.3 
GCD-CPNT(I) 

FHOPS  -  IFHOPS 
F2MUF  -  SF2MUF 


CALL  ELAYER 
CALL  FLAYER 
CALL  XSLOSS 


k - 1 


THE  COMMON  VARIABLES  ARE: 

/PWR/  PWRWAT 

/MUF2/  SF2MUF.  IFHOPS 

/HITE/  HPF2 

/SYSDAT/  MISC  (64) 

DATA  INITIAUZATION  IS  REQUIRED  FOR: 

R2.  AMINO.  XLZ.  DIST 

INITIAUZATION  IS  REQUIRED  EACH  CALL  FOR: 

_  FOEMAX,  FOEMIN,  FLDST.  TLOSS.  ABPS. 

LOCAL  VARIABLES  THAT  GET  VALUES  FROM  GLOBAL 
VARIABLES  ARE  DESCRIBED  IN  THIS  PROCESS  BLOCK 


FRAT  -  FREQ/EMUF 

ADX  -1.3S9  +8.617  •  LOG{FRAT) 

ADX  -  MIN  (ADX.  0) 


TLOSS(I)  -  TLOSS(I)  +  ADX 

FLDST(1 )  -  107.2  +  PWRDB  +  8.686  *  LOG(FREQ)  -  TLOSS(I) 


CALL  MODLOS 

FLDST(2)  «  107.2  +  PWRDB  +  8.686  *  LOG(FREQ)  -  TLOSS(2) 


FHOPS -  0 


FHOPS  -  FHOPS  +  1 
CALL  GEOHOP 


Figure  2.  Subroutine  LTLFLD  (Sheet  1  of  2) 
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I  IE -NF  MODE 


FRAT  -  FREQ/(FOEMAX  ’SECDEF) 
ADX  -1.359  -t-S.ei  7  *  LOG  (FRAT) 
ADX  -  MIN  (ADX.  0) 


TLOSS(3)  -  TLOSS(3)  +  ADX 

FLDST(3)  -  107.2  +  PWRDB  +  8.686  •  LOG  (FREQ)  -  TLOSS(3) 


DBU  -  MAX(  FLDST(1).FLDST(2).FLDST(3) ) 
OMUF-MAX{EMUF.F2MUF) 


RETURN 


Figure  2.  Subroutine  LTLFLD  (Sheet  2  of  2) 


MUF85.  Return  the  greater  of  the  two  in  variable  OMUF.  Note 
that  in  this  case,  OMUF  is  not  the  "operational  MUF"  as  used  in 
Damboldt's  work,  but  the  higher  of  the  E-  and  F-LAYER  MUFs 
resulting  from  different  prediction  methods. 

Important  variables  local  to  LTLFLD  are  listed  below: 

PWRDB  transmitter  power  (dBW) 

GCD  great  circle  distance  (radians) 

FHOPS  number  of  hops  for  F  mode  propagation 

F2MUF  F  layer  MUF  from  MUF85 

ADX  correction  for  E  layer  penetration 

TLOSS  transmission  loss 

FLDST  predicted  field  strength 

HQPEF  number  of  hops  for  mixed  E-F  propagation 

HTEF  weighted  average  height  mixed  mode 

EMUF  E-layer  MUF  as  predicted  in  ELAYER 

ABSI  This  variable  occurs  in  several  subroutines,  with  the 

same  name,  and  is  the  absorption  index 
FQEMIN  minimum  FOE  found  along  path 

FOEMAX  maximum  FOE  found  along  path 

SECD  the  secant  of  the  angle  of  incidence  of  a  ray  path  with 

D  region  ionization  height  (100  km)  for  E  and  F  modes 
FHOPS  the  number  of  F  mode  hops 

SECPHE  secant  of  phi  with  the  E,F  layers 

SECDEF  secant  of  phi  with  the  D  region  ionization  for  mixed  modes 

2.2.4. 1.3  LTLFLD  output s 

DBU  field  strength  (dBU>luV/M) 

OMUF  classical  MUF  (MHz)  (E  or  F2-layer) 

2. 2. 4. 2  Subroutine  ELAYER 

2. 2. 4. 2.1  ELAYER  inputs 
Call  variables 

SSN  sunspot  number  (Zurich  running  average) 

SCD  great  circle  distance,  transmitter  to  receiver  (KM) 

FOEMIN  minimum  FOE  of  all  control  points  used  (MHz) 

FOEMAX  maximum  FOE  of  all  control  points  used  (MHz) 

Common  variables 

CONNUM  the  number  of  control  points  in  use 

D2R  factor  for  converting  degrees  to  radians 

2. 2. 4. 2. 2  ELAYER  processing 

Figure  3  illustrates  the  flow  of  subroutine  ELAYER.  The 
primary  purpose  of  this  subroutine  is  to  compute  the  values  of 
the  E-layer  MUF  (EMUF).  Other  associated  outputs  are  listed  in 

2  4  2 

^  m  ^  ^  m  Wa 
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ELAYER 


DECLARATIONS 

AND 

INITIALIZATION 


SUMABS -  0 


[the  common  variables  used  ARE; 

_ I  /POINTS/CONNUM 

I  /PIDEGS/D2R 

[^TA  INITIALIZATION  REQUIRED  FOR  RZ.  AFLAG 


EKM.1 


EKM>CONN^ 


EKM  -  EKM  +  1 


IKM-EKM 

CALL  GETDAT 

CHI-MIN(CHI.102*D2R) 

ABSUAINDEX(CHI,SSN) 

E2000-[I[[[-4.99*ABSI+29.5irABSI-70.97]*ABSI+89.97TABSI-68.07] 

•ABSI+38.43]*ABSI+3.41 

FOE-E2000*.2085  (E-Layer  distance  factor  for  distance  •  0) 


FOE>FOEMAX 


FOEMAX  -  FOE 


FOE<FOEMIN 


FOEMIN-FOE 


SUMABS  -  SUMABS  -t-  ABSI 
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Rgure  3.  Subroutine  ELAYER  (Sheet  2  of  2) 
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Significant  local  variables  are  listed  below: 

CHI  solar  zenith  angle 

AESP  An  average  value  for  ABSI  across  the  path 

FQE  E  layer  critical  frequency 


2. 2. 4. 2. 3  ELAYER  outputs 

EHOPS  number  of  E-mode  hops  for  this  geometry 

ABSI  absorption  index 

FOEMIN  minimum  FOE  of  all  control  points  used  CMHz) 

FOEMAX  maximum  FQE  of  all  control  points  used  (MHz) 

EMUF  E- layer  MUF  CMHz) 

DELE  take-off  angle  associated  with  this  mode  (deg) 

SECE  secant (phi )  at  E-layer  (110  KM)  for  this  DELE 

DGYRO  9y>'o  frequency  at  D  height  (100)  (MHz) 

SECD  secant (phi)  at  D-layer  (100  KM)  for  this  geometry 

(KM) 

PATHE  total  slant  distance  of  the  path  for  this 

geometry  (KM) 

ABSP  average  absorption  index 

2. 2. 4. 3  Subroutine  FLAYER 

2. 2. 4. 3.1  FLAYER  inputs 
Call  variables 

FHOPS  number  of  F-mode  hops  for  this  geometry 

GCD  great  circle  distance,  transmitter  to  receiver  (KM) 

HPF2  virtual  height  of  F2-layer  (KM) 

F2MUF  F2-layer  MUF  (MHz) 

SSN  sunspot  number  (Zurich  running  average) 

FREQ  current  frequency  (MHz) 

MONTH  current  month 

Common  var i ab 1 es 

CONNUM  the  number  of  control  points  in  use 

CONLAT  array  of  control  point  geographic  latitudes  (deg  N  &  S) 

CONLMT  array  of  control  point  local  mean  times  (hours:  0:23) 

2. 2. 4. 3. 2  FLAYER  processing 

Figure  4  illustrates  the  flow  of  subroutine  FLAYER.  The 
primary  purpose  of  this  subroutine  is  to  compute  the  value  of  the 
critical  frequency  of  the  F2-layer,  F0F2.  Other  associated 
ouptuts  are  listed  in  2. 2. 4. 3. 3. 

There  are  no  significant  local  variables. 
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FLAYER 


DECLARATIONS 

AND 

INITIALIZATION 


CALL  GEOHOP 
FOF2  -  F2MUF/SECF 
INDCON  .  (CONNUM/2)  +  .5 
CALL  F2DIS 


RETURN  ^ 


THE  COMMON  VARIABLES  USED  ARE: 
/POINTS/CONNUM 


Figure  4.  Subroutine  FLAYER 
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2. 2. 4. 3. 3  FLAYER,  outputs 


SECF  secant (phi)  at  F2-height  for  this  geometry 

F0F2  critical  frequency  of  F2-layer  (MHz) 

SIG  standard  deviation  of  F2-layer  MUFS 

SECD  secant (phi)  at  D  height  (100  KM) 

DELF  take-off  angle  for  this  geometry  (deg) 

PSI  included  angle  from  transmitter,  earth  center, 

receiver  (deg) 

PATHF  total  slant  range  of  path  for  this  geometry  (KM) 

2. 2. 4. 4  Subroutine  XSLOSS 
2. 2. 4. 4.1  XSLOSS  inputs 
Call  variables 

GCD  great  circle  distance,  transmitter  to  receiver  (KM) 

Common  variables 


CONNUM  number  of  control  points  in  use 

2. 2. 4. 4. 2  XSLOSS  processing 

Figure  5  illustrates  the  flow  of  subroutine  XSLOSS.  This 
routine  compute  the  excess  loss  adjustment,  ADJ. 

Significant  local  variables  are  listed  below: 

ASM  median  excess  system  loss 

2. 2. 4. 4. 3  XLOSS  outputs 

ADJ  excess  loss  adjustment 

2. 2. 4. 5  Subroutine  MODLOS 
2. 2. 4. 5.1  MODLOS  inputs 
Call  variables 


FREQ 

DGYRO 

ABSI 

HOPS 

PATH 

SECD 

DIST 

LAYER 

MUF 

FO 

SECPHE 

ADJ 


current  frequency  (MHz) 

gyro  frequency  at  100  KM  (MHz) 

absorption  index 

number  of  hops  in  this  mode 

slant  range  of  this  mode  (KM) 

secant (phi)  at  D  height 

distribution  constant 

layer  index,  1=E,  2=F 

maximum  usable  frequency  this  layer  (MHz) 
critical  frequency  (MHz)  -  not  currently  referenced 
secant (phi)  -  not  currently  referenced 
excess  system  loss  adjustment  from  XSLOSS  (dB) 
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Figure  5.  Subroutine  XSLOSS 
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2. 2. 4. 5. 2  MODLQS  processing 

Figure  6  illustrates  the  flow  of  subroutine  MODLOS.  This 
routine  computes  transmission  losses,  TLOSS. 

Significant  local  variables  are  listed  below: 

GZ  Ground  loss 

FSLOSS  Free  space  loss 

ABPS  deviative  absorption 

XLS  additional  loss  due  to  propagation  at  frequency  >  MUF 

2. 2. 4. 5. 3  MODLQS  outputs 

TLOSS  transmission  losses  CdB) 

2. 2. 4. 6  Subroutine  GEOHOP 

2. 2. 4. 6.1  GEOHQP  inputs 
Call  variables 

HOPS  the  number  of  hops  for  which  geometry  is 

determi ned 

GCD  great  circle  distance  (radians) 

HEIGHT  height  of  layer  for  this  mode  (KM) 

2. 2. 4. 6. 2  GEOHOP  processing 

Figure  7  illustrates  the  flow  of  subroutine  GEOHOP.  The 
purpose  of  this  routine  is  to  compute  the  values  of  the  angles 
and  distances  associated  with  the  geometry  of  a  hop. 

2. 2. 4. 6. 3  GEOHQP  outputs 

SECPHE  secant (phi)  at  current  layer  height 

DEL  take-off  angle  for  this  geometry  (deg) 

PATH  total  slant  range  for  this  geometry  (KM) 

PSI  included  angle  for  half  hop,  radians  from  earth 

center  (deg) 

SECD  secant (phi)  at  D-layer  (100  KM)  height  (deg) 

2. 2. 4. 7  Function  AINDEX 
2. 2. 4. 7.1  AINDEX  inputs 

Call  variables 

CHI  solar  zenith  angle  (deg) 

SSN  sunspot  number  (Zurich  running  average) 
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MODLOS 


Figure  6.  Subroutine  MODLOS 
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GEOHOP 


I  DATA  INITIALIZATION  IS  REQUIRED  FOR 
1  R2.  RI02 


PSI  -  .5  *GCD/HOPS 
COSPSI  -  COS  (PSI) 

SINPSI .  SIN(PSI) 

TANPHE  -  SINPSI/(1 -COSPSI  +  HEIGHT/RZ) 

PHE  -  ATANfTANPHE) 

SECPHE-  1/COS{PHE) 

DEL  -  PI02-PHE-PSI 
COSDEL  -  COS(DEL) 

SINDEL  -  SIN(DEL) 

SINPHED  -  R2  •COSDEL/(R2  +110) 

SECD  -  1/SQRT(1-S1NPHED**2) 

PATH  .  ABS(HOPS  *  2  *  (HEIGHT  +  R2‘  ,1-COSPS''1*SECPHE) 


RETURN  ^ 


DECLARATIONS 

AND 

INITIALIZATION 


Figure  7.  Subroutine  GEOHOP 
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2. 2. 4. 7. 2  AINDEX  processing 

Figure  8  illustrates  the  flow  of  function  AINDEX. 
function  computes  the  value  of  the  absorption  index,  AINDEX. 

2. 2. 4. 7. 2  AINDEX  outputs 

AINDEX  absorption  index 

2. 2. 4.8  Subroutine  SYSSY 

2. 2. 4. 8.1  SYSSY  inputs 
Call  variables 

SLAT  geomagnetic  latitude  of  reflection  area 

(deg  N  S> 

Z  reflection  area  local  mean  time  (hours:  0-23) 

GCDRAD  great  circle  distance  (radians) 

Common  variables 

R2D  factor  for  converting  radians  to  degrees 

2. 2. 4. 8. 2  SYSSY  processing 

Figure  9  illustrates  the  flow  of  subroutine  SYSSY.  This 

routine  computes  the  excess  system  loss,  FM. 


2. 2. 4. 8. 3  SYSSY  outputs 

FM  excess  system  loss  (dB) 

2. 2. 4. 9  Subroutine  F2DIS 

2. 2. 4. 9. 1  F2DIS  inputs 
Call  variables 

FMUF  F2MUF  (MHz) 

SSN  sunspot  number  (Zurich  running  average) 

CLAT  geographic  latitude  (deg  N  &  S) 

FREQ  operating  frequency  (MHz) 

AB  local  mean  time  (hours:  0-23) 

MONTH  current  month 

2. 2. 4. 9. 2  F2niS  processing 

Figure  10  illustrates  the  flow  of  subroutine  F2DIS.  This 
routine  computes  the  standard  deviation  of  the  normal  F2MUF. 
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AINl 

_ ^ 

3EX  ^ - [  REAL-TYPED  FUNCTION 

ABSl  -  (1  +  .0037  ‘SSN)  *{C0S(.881  *CHI))  **1.3 

ABSIMN  .  .025  *  (1  -t-  .031  *SSN) 
AINOEX  -  MAX  (ABSl.  ABSIMN) 

_j  CHECK  ABSl  AGAINST  MINIMUM 

l^i^wed  nighttime  absorption 

URN  ^ 

Figure  8.  Function  AINDEX 
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Figure  1 0.  Subroutine  F2DIS  (Sheet  1  of  4) 
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Figure  10.  Subroutine  F2DIS  (Sheet  2  of  4) 
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E 


Figure  10.  Subroutine  F2DIS  (Sheet  3  of  4) 
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Figure  1 0.  Subroutine  F2DIS  (Sheet  4  of  4) 
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CALXLS 


Figure  1 1 .  Function  CALXLS 
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2.2.4.11  Subroutine  GETDAT 


2.2.4.11.1  GETDAT  inputs 
Call  variables 

IKM  control  point  index  (l=first,  2=second,  etc.) 

Common  variables 


CONLAT 

array  of  control 
(deg  N  &  S) 

poi  nt 

geographic 

1  at i tudes 

CONLON 

array  of  control 
(deg  E  &  U) 

poi  nt 

geographic 

longitudes 

CONLMT 

array  of  control 
(hours:  0:23) 

point 

local  mean 

times 

CONMLAT 

geomagnetic  latitudes 
(deg  N  &  S) 

of  control 

points  in  CPMT 

2.2.4.11.2  GETDAT  processing 

Figure  12  illustrates  the  flow  of  subroutine  GETDAT.  This  is 
routine  translates  variable  names  and  formats  of  data  for 
compatibility  with  the  rest  of  the  field  stength  module. 

2.2.4.11.3  GETDAT  outputs 

LO  geographic  longitude  of  control  point  IKM  (deg  E  &  Ui) 

UO  geographic  latitude  of  control  point  IKM  (deg  N  &  S> 

LMTNOW  local  mean  time  of  control  point  IKM  (hours:  0-23) 

MLAT  geomagnetic  latitude  of  control  point  IKM  (deg  N  &  S) 

CHI  solar  zenith  angle  of  control  point  IKM  (deg) 

2.2.4.12  Function  ZENANG 

2.2.4.12.1  ZENANG  inputs 
Call  variables 

LO  geographic  longitude  (deg  E  &  W) 

WO  geographic  latitude  (deg  N  &  S) 

COMMON  var i ab 1 es 

SLAT  subsolar  latitude  (deg  N  &  S) 

SLON  subsolar  longitude  (deg  E  &  W) 

2.2.4.12.2  ZENANG  processing 

Figure  13  illustrates  the  flow  of  function  ZENANG.  This 
function  computes  the  zenith  angle,  ZENANG. 
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GETDAT 


r 


I 


L 


LO  -  CONLAT  (IKM) 

WO  -  CONLON  (IKM) 
LMTNOW  .  CONLMT  (IKM) 
MLAT  -  CONMUVT  (IKM) 
CHI  .  ZENANG  (LO.  WO) 


RETURN^ 


DECLARATIONS 

AND 

INITIALIZATION 


THE  COMMON  VARIABLES  USED  ARE: 

/POINTS/CON  LAT,  CONLON.  CONLMT,  CONMLAT 
NO  DATA  INITIALIZATION  IS  REQUIRED 


Figure  12.  Subroutine  GETDAT 
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ZENANG 


I  THE  COMMON  VARIABLES 

- 1  USED  ARE: 

I  /SUN/SUT,  SLON 

(mw  m 


CHI  -  SIN{L0)  •SIN(SLAT)  +  COS(LO)* 
COS(SLAT)  •COS{WO-SLON) 
CHI -MAX  (CHI,-1.) 

CHI  -  MIN  (CHI.  1) 

ZENANG  -  ACOS  (CHI) 


RETURN^ 


DECLARATIONS 

AND 

INITIALIZATION 


Figure  13.  Function  ZENANG 
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2.2.4.12.3  ZENANG  outputs 

ZENANG  zenith  angle  (deg  N  &  S) 

2.2.4.13  Function  FNORMAL 

2.2.4.13.1  FNORMAL  inputs 

Y  number  of  standard  deviations  from  distribution  center 

2.2.4.13.2  FNORMAL  processing 

Figure  14  illustrates  the  flow  of  function  FNORMAL.  This 
function  computes  the  cumulative  normal  distribution,  FNORMAL, 
from  a  maximum  of  5  standard  deviations. 

Significant  local  variables  are  described  below: 

Y  number  of  srtandard  deviations  above  the  mean 

S  cumulative  normal  distribution  from  the  mean 

2.2.4.13.3  FNORMAL  outputs 

Y  number  of  standard  deviations  from  distribution 

center 

FNORMAL  cumulative  normal  distribution  (from  maximum  of 
5  standard  deviations) 
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Propaqation  Enqineerinq  Division,  Technical  Report,  EMEO-PED-79- 
7,  Sep  1978. 
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Including  an  Ionospheric  Model,  NRL  Report  222L,  Mar  1971. 
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FNORML 


DECLARATIONS 

AND 

INITIALIZATION 


_ j_  NO  COMMON  VARIABLES  ARE  USED. 

I  DATA  INITIALIZATION  IS  REQUIRED  FOR  SQRT2 

L_. 


Figure  14.  Function  FNORML 
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9.  K.  A.  Norton,  Transmission  Loss  in  Radio  Propagation,  NBS, 
Tech  Report  Note  12,  Jun  1959. 


2.3  SOUNDER  UPDATE  (SU) 

2.3.1  Sc  ope 

2.3. 1.1  Identification 
SUBROUTINE  SOUNDR 

Callable  through  the  existing  module  REal-time  (GETRTM)  by 
first  entering  command  "RE“  through  the  main  line  "TDA"  program, 
and  then  responding  'Y’  to  the  question  "DO  YOU  WANT  REAL  TIME 
DATA";  or  by  the  new,  direct  command  "SU"  enter able  through  the 
main  line  "TDA"  program. 

2.3. 1.2  Programmer 

D.  J.  Brandon  —  Systems  Exploration,  Inc. 

2.3. 1.3  Purpose 

The  sounder  update  module  allows  the  user  to  calculate  a 
’local'  pseudo  Sun  Spot  Number  (SSN)  and  pseudo  flux  value,  based 
on  the  location  of  a  user-defined  transmitting  and  receiving 
antenna,  and  a  reference  sounder  path  MOF,  which  is  then  used  to 
derive  a  more  accurate  MUF  prediction  for  other  paths. 

2.3. 1.4  Background 

The  computer  model  development  for  Sounder  Updates  in  this 
system  is  a  re-hosted  version  of  the  BASIC  model  operating  on  the 
Tektronix  4052  Graphics  System  for  the  Army  Prophet  Evaluation 
System  (APES). 

Reference  to  external  routines  were  modified  to  match  those 
used  in  the  TDA  operating  environment,  where  a  similiar  function 
already  existed.  New  subroutine  CVTLL  was  written  to  handle  the 
numeric  conversions  required  to  list  out  the  site  latitudes  and 
longitudes  as  that  function  was  not  readily  available  in  the 
existing  TDA  subroutine  library. 
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^  m  ^  9  X  • 


Restrictions  and  Limitations 


The  sounder  update  module  is  built  around  the  MUF85  model, 
with  a  MDF  range  between  2  and  50  which  is  consistent  with  other 
models  operating  within  the  TDA  environment. 

The  new  model  was  enhanced  to  solve  for  the  computed  MUF 
value  to  three  decimal  points,  increase  SSN  range  up  to  a  value 
of  250,  and  be  able  to  converge  on  the  proper  MUF  values  from 
either  direction. 

2.3. 1.6  Language 
HP  FORTRAN/77 

2.3. 1.7  Computer  Configuration 

HP  SOOO  Series  500  HP-UX  (HP  9050) 

2.3. 1.8  Models  Referenced 

AS,  DA,  RE,  SS,  TDA,  TM 

2.3.2  Numeric  Method 

The  SORFMS  program  will  be  able  to  use  information  obtained 
from  an  HF  oblique  sounder  to  update  the  value  of  10.7  cm  flux 
used  in  the  MINIMUF  model.  To  take  advantage  of  this  feature,  a 
sounder  path  must  be  defined  and  the  maximum  observed  freqeuncy 
(MOF)  from  the  sounder  must  be  updated  periodically.  The  user 
enters  receiver  and  transmitter  parameters  associated  with  the 
sounder  site  including  latitude,  longitude,  antenna,  power,  and 
antenna  orientation.  Once  these  parameters  are  established  for  a 
given  sounder  path,  the  user  may  skip  this  step.  The  user  enters 
the  current  MOF  for  the  sounder  path.  The  system  then  will 
compute  a  value  of  10.7-cm  flux  that  will  cause  the  computed  MUF 
for  the  sounder  path  to  agree  with  the  MOF.  This  value  of  10.7- 
cm  flux  will  then  be  used  for  all  MUF  calculations  until  it  is 
changed  by  another  SU  command,  or  until  a  different  value  is 
specified  by  the  FL  (flux)  command.  To  provide  this  sounder 
update  capability  in  SORFMS,  the  Sounder  Update  module  was 
rehosted  from  APES. 

2.3.3  Module  Design  Description 
2.3.3. 1  Required  Input 

a.  User  entries  to  initiate  the  Sounder  Update  model : 

1.  Date  of  the  sounder  MOF  specified  stored  in  variables 
MONTH,  DAY,  JDAY,  YEAR.  (Set  by  the  TDA  command  'DA».) 

2.  Time  of  the  sounder  MOF  specified  stored  in  variables 
HOUR  and  MINUTE.  (Set  by  the  TDA  command  'TI'.) 
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3.  Sun  Spot  Number  (SSN>  of  the  sounder  MOF  specified  stored 
in  variable  (TBS).  (Set  by  the  TDA  command  ’SS’.) 

4.  Transmitter  site,  including  the  site  latitude 
(STALAT (KXMT) )  and  longitude  (STALON (KXMT) ) .  (These  values  may  be 
set  up  by  the  TDA  command  'AS^,  but  only  the  site  Name  may  be 
entered  within  the  Sounder  module.) 

5.  Receiver  site,  including  the  site  latitude  (STALAT (KRCV) ) 
and  longitude  (STALON (KRCV) .  (These  values  may  be  set  up  by  the 
TDA  command  ’AS’,  but  only  the  site  Name  may  be  entered  within 
the  Sounder  module.  ) 

6.  Sounder  MOF  on  this  path  (CMUF). 
b.  Subroutine  SOUNDR  Inputs 

No  called  variables 
COMMON  Variables: 


LUNI  CRT  input  logical  unit  (5) 

LUNO  CRT  output  logical  unit  (6) 

MNEM  Array  of  three— char acter  months  of  the  year 

DTR  Conversion  factor  for  degr ees-to-r adi ans 

LISTAl  Primary  array  of  input  characters  received  from  user 

LISTBl  Secondary  array  of  input  characters  received  from  user 

MONTH  Month  of  model  run 

DAY  Day  of  model  run 

MINUTE  Minute  of  model  run 

NUMSTA  Number  of  stations  (sites)  available  in  TDA  environment 
STAANT  Array  of  antenna  number  (type)  for  all  sites 
FLXIO  10.7  cm  flux 

SSN  Sun  spot  number 

STALAT  Array  of  latitudes  for  all  sites 

STALON  Array  of  longitudes  for  all  sites 
STABRG  Array  of  bearings  for  all  sites 
STAPWR  Array  of  powers  for  all  sites 

2. 3. 3. 2  Processing 


Figures  15  and  IS  illustrate  the  flow  of  the  SU  module.  Once 
the  model  system  is  initialized,  the  user  is  prompted  to  enter 
the  transmitter  and  receiver  sites  to  define  the  path  in 
question.  He  is  then  prompted  to  enter  the  sounder  MOF  of  the 
specified  path.  These  values,  along  with  the  system  date,  time 
and  SSN  are  then  used  to  convert  the  supplied  MOF  to  the  MUF 
value  of  this  instant. 


THE  COMMON  VARIABLES  USED  ARE 


(^SOUNDER ^ 


■1 

DECLAR 

AN 

INITIALI2 

ATIONS 

D 

ATION 

/DATIME/MNEM 
/PIDEGS/DTR 
/STRING/LISTA.  LISTS 
/SYSDAT/STANAM,  MONTH.  DAY,  YEAR. 

JDAY.  HOUR.  MINUTE.  NUMSTA. 

I  STAANT,  MISC  (1-6,  STABRG.  STAPWR 

■f  FLX1 0,  SSN,  STALAT,  STALON 

I  /OPRDEV/  LUNI,  LUNO,  PROMPT 
I  /HITE/V 

I  DATA  INITIALIZATION  IS  REQUIRED  FOR 
I  SMOFLO,  SMOFHI,  SSNMAX 

I  INITIALIZATION  IS  REQUIRED  EACH  CALL  FOR 
I  PROMPT.  RESU.  V 


Figure  15.  Subroutine  SOUNDR  (Sheet  1  of  9) 
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0 


Figure  15.  Subroutine  SOUNDR  (Sheet  2  of  9) 
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Figure  15.  Subroutine  SOUNDR  (Sheet  3  of  9) 
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YES 


DISPLAY  MOF  ON  CRT  (SMOF) 

CALL  CVTLL  FOR  TRANSMITTER  LATITUDE  AND  LONGITUDE 

DISPLAY  SITE  NAME.  LAT.  LONG.  ANTENNA  #.  POWER.  BEARING  FOR  TRANSMITTER 

CALL  CVTLL  FOR  RECEIVE!"  LATITUDE  AND  LONGITUDE 

DISPLAY  SITE  NAME.  LAT  LONG.  ANTENNA  #.  POWER.  BEARING  FOR  RECEIVER 

CONVERT  TRANSMITTER  LATITUDE  TO  RADIANS  (TLAT) 

CONVERT  TRANSMITTER  LONGITUDE  TO  RADIANS  (TLON) 

CONVERT  RECEIVER  LATITUDE  TO  RADIANS  (RLAT) 

CONVERT  RECEIVER  LONGITUDE  TO  RADIANS  (RLON) 


Figure  15.  Subroutine  SOUNDR  {Sheet  4  of  9) 


Figure  15.  Subroutine  SOUNDR  (Sheet  5  of  9) 
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SMOF  LO  <CCMUF<SMOF  HI 
AND 

CMUF  •GE’  SMOF 
s.  AND  ^ 

X.SSN  IE’  SSN  MAX/ 


I  IS  CALCULATED  MUF  WITH  IN 
■  j  RANGE  AND  >  OBSERVED  MOF 


SSN.SSN.»-10 
CALL  MUF85 


^  IS  ^ 
CMUF  <  SMOF 


DIRFLG-  -1 


IS  CALCULATED  MUF 
<  OBSERVED  MOF 


SSN-SSN  +  (DIRFLG  *  STEPSZ)  /2 


CALL  MUF85 


'  ADJUST  SSN  TO  MIDPOINT 
j  OF  SOLUTION  INTERVAL 


Figure  15.  Subroutine  SOUNDR  (Sheets  of  9) 


Figure  15.  Subroutine  SOUNDR  (Sheet  7  of  9) 


Figure  15.  Subroutine  SOUNDR  (Sheet  8  of  9) 


A7 


YES  ^ 

NFMT  .  7 

YES 


NFMT-NELEM 


I  DISPLAY  A  LIST  OF  THE 
-1  AVAILABLE  SOUNDER  SITES 
I  UP  TO  7  SITES  PER  LINE 
I _ 


RETURN  TO  OPERATOR 
■  PROMPT  WHERE 
WAS  ENTERED 


NELEM  -  NELEM-7 
INCB-r  (M)  +1 
INCE-r  (1-1)  +NFMT 


WRITE  (LUNO,  7) 


{  DISPLAY  UP  TO  7 
SITE  NAMES 


Figure  15.  Subroutine  SOUNDR  (Sheet  9  of  9) 


CONVERT  LATITUDE  AND  LONGITUDE 
TO  A  PRINTABLE  FORMAT 
CONSISTANT  WITH  APES 

LATITUDE  PRINTS  ‘S’  IF  NEGATIVE 
•N’  IF  POSITIVE 

LONGITUDE  PRINTS  ’E’  IF  NEGATIVE 
-  W  IF  POSITIVE 


Figure  16.  Subroutine  CVTLL 


The  Sun  Spot  Number  CSSN)  and  related  10.7  cm  flux 
values  are  obtained  by  adjusting  the  computed  MUF  for  the  sounder 
path,  by  varying  the  SSN  until  it  agrees  with  the  supplied  MOF, 
or  until  50  iterations  have  been  completed. 

The  SSN  is  then  reset  to  zero  and  then  cont*  .^aiiy  increased 
by  units  of  10  until  the  resultant  MUF  value  realized  by  the 
MUFS5  model  calculations  exceed  the  supplied  MOF  value,  or  until 
the  SSN  exceeds  250. 

The  SSN  is  further  adjusted,  by  performing  a  binary 
iteration  on  the  resultant  MUF  value  obtained  from  the  MUF85 
model,  until  the  difference  between  this  MUF  value  and  the 
originally  supplied  MOF  value  is  .0005,  or  until  50  iterations 
have  occured.  The  10.7  cm  flux  value  is  then  computed,  based  on 
this  new  SSN. 

The  original  and  final  SSN’s  and  10.7  cm  flux  values  are 
printed  out.  The  user  is  prompted  as  to  whether  or  not  the  final 
values  should  be  used  to  replace  the  orginal  values  in  permanent 
TDA  memory,  to  be  used  by  this  and  all  other  TDA  modules. 

The  user  is  then  prompted  to  repeat  these  calculations  using 
a  different  MOF  at  the  same  sites.  If  he  responds  with  a  ’y' 
then  the  entire  process  repeats;  otherwise  the  sounder  update 
module  terminates,  and  control  is  returned  to  the  calling  process 
(either  the  REal-time  module  or  the  TDA  main  program) 

2. 3. 3. 3  Output 

1.  The  updated  SSN  value  (if  the  user  responds  'y'  to  the 

update  request) 

2.  The  corresponding  10.7  cm  flux  value 


2.3.4  References 

1.  NOSC  TD  782,  Tactical  Decision  Aids  for  HF  Communication,  by 
D.B.  Dailors,  15  November  1S84. 

2.  NOSC  TD  10&5,  Operational  User's  Manual  for  Tactical  Decision 
Aids  for  HF  Communication,  by  D.B.  Sailors,  July  1987. 

3.  NOSC  Contractor  Report  CR  292,  Operational  Users  Manual  for 
Army  PROPHET  Evaluation  System  (APES) /Theater  Nuclear  Forces 
(TNF)  Frequency  Managment  System,  prepared  by  J.R.  Gnessin, 
Systems  Exploration,  Inc.,  July  1985. 

4.  NOSC  Contractor  Report  (unpublished).  Program  Package 
Document  for  the  Army  PROPHET  Evaluation  System  (APES) /Theater 
Nuclear  Forces  (TNF)  Freqeuncy  Management  System,  prepared  by 
J.R.  Gnessin,  Systems  Exploration,  Inc.,  Oct  1985. 
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5.  HP  FORTRAN/77  Reference  Manual,  Hewlett-Packard  97081-90010, 
June  1985 

.  HP-UX  Reference,  HP9000  Series  500  Computers,  HP-UX  Release 
.2,  Hewlett-Packard  09000-90010,  Edition  1,  April  1987. 


3.1  SU  TEST 


SECTION  3 
TESTING 


3.1.1  Purpose 

The  purpose  of  testing  the  SU  (nodule  is  to  verify  that  the 
SU  module,  as  translated  from  the  APES  Tektronix  BASIC  to  TDA  HP 
9050  FORTRAN,  operates  in  essentially  the  same  manner  in  both 
systems.  Because  of  differences  in  the  method  of  MUF 
computations  in  the  two  systems,  intermediate  results  cannot  be 
effectively  compared  —  namely  the  pseudo  sunspot  numbers 
computed  under  the  different  methods.  However,  the  purpose  of 
the  SU  module  is  to  shift  the  MUF  curve  so  that  the  MUF  equals 
(or  approximately  equals)  the  MOF  at  the  measured  time.  Thus, 
the  purpose  of  the  test  is  to  demonstrate  that  SU  module  performs 
this  function  in  the  HP  9050  TDA  and  differences  in  the  pseudo 
sunspot  numbers  are  considered  irrelevant. 

J.1.2  Procedure 

Testing  the  SU  module  consisted  of  entering  selected  MOF's 
and  verifying  that  the  SU  module  operates  as  intended.  The  data 
was  for  a  sounder  path  from  Isabela,  Puerto  Rico  to  Norfolk,  VA 
on  15  and  IS  November,  1981.  Data  was  obtained  from  the 
following  report:  NRL  Memorandum  Report  5284,  H.F.  Frequency 
Management  by  Frequency  Sharing  as  Assisted  by  Models  Updated  in 
Real-Time,  by  D.R.  Uffelman  and  L.O  Harnish,  SRI  International, 
and  J.M.  Goodman,  E.Q.  Hulburt  Center  for  Space  Research, 
Ionospheric  Effects  Branch,  Space  Science  Division,  March  21, 
1984. 


3.1.3  Results  and  Conclusions. 

In  general,  when  a  MOF  value  was  entered  into  the  system  via 
the  SU  module  prompt,  the  SU  module  produced  a  pseudo  sunspot 
number.  Then,  when  the  pseudo  sunspot  number  was  entered  into 
the  system  to  replace  the  previous  sunspot  number,  the  system 
computed  a  corresponding  MUF  which  was  equal  or  nearly  equal  to 
the  original  MOF. 

The  test  demonstrated  that  the  MOF’s  on  the  particular  days 
and  paths  selected  stressed  the  model  to  its  upper  limits.  The 
model  provides  for  sunspot  numbers  up  to  250.  In  many  cases, 
using  measured  MOF’s  for  the  selected  path  and  times,  the 
computed  pseudo  sunspot  number  exceeded  the  upper  limit.  For 
those  MOF’s  which  were  in  range  (approximately  15  MHz  or  less 
during  early  day  time  hours)  the  model  worked  as  intended.  When 
the  MOF  was  out  of  range,  the  model  could  not  produced  an 
equivalent  MUF. 

The  bounding  contraints  of  the  SU  model  are  included  by 
design  to  avoid  ambiguous  results.  The  physical  relationship  of 
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MUF  to  sunspot  number  is  that  MUF  increases  with  increasing 
sunspot  number  to  a  maximum  value  (saturation)  and  then  begins  to 
Thus,  the  model  cuts  off  further  computations  at  the 
saturation  level.  The  SU  code  was  modified  to  test  for  HOF’s 
above  the  maximum  allowable  value  for  the  given  path  and  time  and 
provide  a  prompt  to  the  user  when  the  value  was  exceeded. 

The  model  performed  as  intended  for  those  MDF’s  which  were 
in  range.  For  comparison,  selected  MOF's  were  also  entered  on 
the  Tektronix  4052  containing  the  APES  model,  and  results  were 
comparable  (though  not  identical  due  to  differences  in  versions 
of  MUF  85  used  in  the  models). 


3.2  FS  TEST 

3.2.1  Purpose 

The  purpose  of  testing  the  FS  module  is  to  verify  that  the 
predicted  field  strength  outputs  of  the  upgraded  model  are  more 
closely  aligned  with  observed  data  than  those  of  the  previous 
version  of  the  model.  This  would  then  demonstrate  that  the 
enhancements  to  the  model  actually  improved  the  accuracy  of  the 
model.  It  is  generally  presumed  that  if  the  accuracy  improved, 
the  physical  assumptions  used  in  the  enhanced  model  were  correct. 
Conversely,  if  the  accuracy  of  the  predicted  outputs  were 
degraded,  then  the  assumptions  would  have  been  incorrect.  This 
also  presumes  that  the  observed  data  was  correctly  recorded. 

The  CCIR  data  base  was  used  to  provide  the  observed  data 
which  was  then  used  for  comparison  with  predicted  outputs 
obtained  from  the  model. 

3.2.2  Procedure 

Figure  17  illustrates  an  overview  of  the  sequence  of  testing 
activities  for  the  FS  module  with  approximate  times  required  to 
accomplish  each  major  task.  Testing  activities  are  described 
below: 

a.  The  FS  module  was  first  developed  and  tested  on  the  VAX. 
Modifications  were  incorporated  as  described  earlier  in  this 
report.  The  upgraded  model  was  then  run  on  the  VAX  and  data 

outputs  were  evaluated  for  general  credibility.  This  was 
accomplished  by  manually  comparing  the  predicted  FS  data  results 
for  long  paths  with  selected  data  available  from  the  previous 
version  of  the  FS  model.  For  short  path  lengths,  no  similar  data 
comparison  was  possible  as  this  was  outside  the  region  of 
operability  of  the  previous  model.  However,  various  manual 

methods  were  used  to  determine  general  credibi 1 itiy,  such  as 
reviewing  published  reports  and  performing  manual  calculations  to 
determine  conformance  with  theory  and  selected  empirical  data. 
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Figure  1 7.  Sequence  of  Activities  to  test  Field  Strength  Module 


b.  The  FS  module  was  then  integrated  into  the  HP  9050  version  of 

TDA  and  run  in  a  standard  operator -inter active  mode.  Visual 
observations  of  the  outputs  on  the  CRT  were  compared  with  results 
of  similar  runs  made  on  the  VAX  to  determine  that  the  model 
performed  satisfactorily  as  intended.  Note  that  this  test  did 
not  validate  the  accuracy  of  results  as  compared  with  real  world 
data;  only  that  the  HP  9050  version  was  functionally  identical  to 
the  VAX  developmental  model.  This  test  was  necessary  to 
determine  that  the  translation  from  VAX  to  9050  was  correct.  It 
also  tests  the  operational  version  of  the  model  characterized  by 
two  dimensional  tabular /graphic  outputs  of  field  strength 

computed  as  functions  of  time  and  frequency. 

c.  The  FS  module  was  then  modified  to  provide  outputs  suitable 

for  comparison  with  the  CCIR  data  base.  This  modification  was 

also  necessary  to  adapt  the  model  to  accept  a  series  of  inputs 

stored  in  a  command  file  which  would  be  used  to  drive  the  model 
into  producing  the  desired  set  of  outputs  for  comparison. 

Table  1  illustrates  one  page  of  this  command  file.  The  left 
column  is  the  actual  command  file  except  for  bracketed 

explanations  included  here  for  clarification.  The  right  column 
illustrates  the  system  prompts  corresponding  to  the  commands. 
The  command  file  and  the  statistical  computation  files  were 
previously  established  by  NQSC  personnel  for  the  purpose  of 
comparing  previous  versions  of  FS  module  and  other  related 
modules. 

The  most  significant  modification  to  the  test  version  of  the 
FS  module  was  to  change  the  output  for  one  dimensional  data: 
field  strength  versus  time  with  frequency  fixed  to  correspond  to 
one  of  the  frequencies  associated  with  the  measured  CCIR  data. 
Table  2  illustrates  a  sample  of  this  output.  Bracketed  column 
heads  were  added  to  the  table  for  clarification.  The  Predicted 
FS  versus  UT  (universal  time)  data  shown  in  Table  2  generally 
corresponds  to  a  single  rc v  of  the  conventional  FS  module  output. 
Conventional  output  is  a  two  dimensional  table  (FS  as  a  function 
frequency  and  time).  A  difference  is  that  the  frequency  on  the 
modified  FS  module  output  table  is  repesented  as  a  real  decimal 
to  tenths  of  a  MHz,  whereas  choices  on  the  conventional 
frequency-indexed  output  are  limited  to  even  integer  values 
(actually  rounded  internally)  between  2  and  30  MHz.  It  is  noted 
here  that  the  CCIR  data  is  actually  constrained  further  in  that 
it  is  not  avialable  for  all  possible  times  but  only  a  subset. 

d.  The  modified  FS  module  was  then  run  using  the  command  file. 
The  command  file  automatically  provides  the  set  of  inputs 
corresponding  to  the  CCIR  database  observation  parameters  and 
generate  the  set  of  data  outputs  to  be  used  for  comparison  with 
the  CCIR  observed  data.  This  run  takes  approximately  60  hours  to 
complete  and  was  performed  either  in  a  background/batch  mode  or 
run  on  weekends  to  avoid  disrupting  ongoing  operations  at  the 
NOSC  HP  9050  computer  center. 
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Table  1.  Sample  Page  of  FS  Command  File  and  Corresponding 
System  Prompts 


COMMANDS/RESPDNSES  SYSTEM  PROMPTS  (Abbreviated) 

2  C 2= ALPHANUMERIC  CRT]  PLEASE  SELECT  [TERMINAL!  TYPE 

NUMBER 

n  WOULD  YOU  LIKE  AN  INTRODUCTION 

<Y/N) 

cs  [change  station!  :> 

hono  [station  name,  CS> 

transmitter  =  Honolulu! 

49.40  [latitude! 

-6.19  [longitude! 

0  [antenna  number! 

1000  [power! 

[blue  force! 

0  [speed! 

cs  [change  station! 
sdiego  [station  name, 
receiver  »  San  Diego! 

51.07 
“7.16 

0 

1000 


0 


da 

[date  analyzed! 

;> 

8 

[month! 

da> 

15 

[day! 

da> 

84 

[year! 

da> 

fr 

[operating  frequency! 

;> 

6. 

1  [minimum! 

fr> 

[maximum! 

fr> 

ssn 

[sun  spot  number! 

;> 

25 

• 

ss> 

fs 

[field  strength! 

:> 

s 

[signal  strength  table! 

[display  next  12  hours  of 

data 

e 

[end! 

da 

9 

15 

84 

fr 

6.1 
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CS> 

CS> 

CS> 

CS> 

CS> 

CS> 


Table  2.  Data  File  tstnew  Sample 


XMTR:  LAT;  49.4  LON;  -6.2  SSN;  25.  8/1984 

RCVR;  LAT;  51.1  LON;  -7.2 


[Universal  Time] 

[Local  Mean  Time] 

[Predicted  FS] 

[Frequency] 
[LUF] 

[MUF] 


1  1. 
2  2. 

3  3. 

4  4. 

5  5.. 

6  £.' 

7  7.. 

8  8.- 

9  9.- 

10  10.. 
11  11. < 
12  12.. 

13  13.. 

14  14.. 

15  15.. 

16  16.. 

17  17.. 

18  18.. 

19  19.. 

20  20.. 
21  21.. 
22  22.. 
23  23.. 


-1.00000 
-3. 00000 
-10.00000 
-6.00000 
1 . 00000 
19.00000 
25.00000 
27.00000 
27.00000 
27.00000 
30.00000 
30.00000 
30.00000 
27.00000 
27.00000 
27.00000 
27.00000 
26.00000 
24.00000 
19.00000 
5.00000 
1 . 00000 
-1.00000 
. 00000 


[Operational  MUF] 


XMTR; 

RCVR; 


LAT; 

LAT; 


49.4  LON; 
51.1  LON; 


-6.2 

-7.2 


6. 10000 
6. 10000 
6. 10000 
6.10000 
6. 10000 
6. 10000 
6. 10000 
6. 10000 
6. 10000 
6. 10000 
6. 10000 
6. 10000 
6.  lOOOO 
6. lOOOO 
6.  lOOOO 
6. lOOOO 
6.  lOOOO 
6. lOOOO 

6. 10000 

6. 10000 
6. lOOOO 
6. lOOOO 
6. lOOOO 
6. lOOOO 

SSN;  15. 


2. 00000 
2. 00000 
2.00000 
2.00000 
2.00000 
2. 00000 
2. 00000 
2.00897 
2.29286 
2.49826 
2. 63299 
2.70455 
2.71260 
2.65728 
2.53916 
2.35160 
2.08837 
2.00000 
2.00000 
2. 00000 
2.00000 
2. 00000 
2.00000 
2. 00000 


1.82359 
3.54893 
2.91485 
3.21894 
4.29077 
4.94040 
5.37735 
5.68215 
5.89435 
6.02582 
6.11955 
6. 15365 
6. 14320 
6.09132 
5.99958 
5.86812 
5.69551 
5.47800 
5.20800 
4.87000 
4.42816 
4.23241 
3.82012 
3.86374 


8.27617 
7.73089 
6.39008 
7.10139 
9.52552 
10.85344 
11.68901 
12.22018 
12.54024 
12.86912 
13.07596 
13. 17727 
13. 18331 
13.01563 
12.76410 
12.43015 
12.01183 
11.57844 
11.03184 
10.33840 
9.42092 
9.04360 
8. 19799 
8.32732 


9/1984 


0 

. 44440 

-7.00000 

1 

1 . 44440 

-9.00000 

2 

2.44440 

-10.00000 

3 

3.44440 

-14.00000 

4 

4.44440 

-3. 00000 

5 

5.44440 

14.00000 

6 

6.44440 

24.00000 

7 

7.44440 

26.00000 

8 

8.44440 

27.00000 

9 

9.44440 

27.00000 

10 

10.44440 

27.00000 

11 

1 1 . 44440 

26.00000 

12 

12.44440 

26.00000 
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6. lOOOO 

6. 10000 

6. 10000 

6. 10000 

6. 10000 

6. 10000 

6.10000 

6. lOOOO 

6. 10000 

6. 10000 
6.10000 
6. 10000 
6. lOOOO 


2. 00000 
2.00000 
2. 00000 
2.00000 
2.00000 
2.00000 
2. 00000 
2.00000 
2. 17407 
2.41428 
2.56677 
2.64034 
2.63825 


3.12210 
3.00113 
2.86469 
2.53986 
3.62395 
4.70086 
5. 19650 
5.50159 
5.69744 
5.81504 
5.86952 
5.86908 
5.81817 


7.27761 
6.89848 
6.49210 
5.67374 
7.97812 
10.37070 
11.48817 
12. 18809 
12.64832 
12.85561 
12.92175 
12.86649 
12.70108 


Table  3.  Data  File  fsdb.hfbc  Sample 


L  *CData  base  code]  C101=HFBC&4  (used  for  LTLFLD  model)] 

\  ^[Number  of  data  base  fields] 

\  [Total  data  base  records] 

\  \  ^S^-«:^Data  base  title] 

101  10  12277  HFBC84  Field  Strength  model 

[Path  ID] 

[Circuit  Length] 

[Season:  l=Winter,  2=5pring,  3=Summer,  4=Fall] 

[Orientation:  0=other,  l=North/South,  2=East/Uest] 
[Smoothed  SSN] 

[Mid  path  latitude] 

[Mid  path  local  time] 

[Freq/MUF  ratio] 

’[Observed  Field  Strength] 

■[Predicted  Field  Strenqth] 


1 

175 

3 

0 

42 

50.2 

6.4 

1.2 

23.0 

32.0 

1 

175 

3 

0 

42 

50.2 

7.4 

1.0 

22.0 

33.0 

1 

175 

3 

0 

42 

50.2 

8.4 

1.0 

22.0 

31.0 

1 

175 

3 

0 

42 

50.2 

9.4 

.9 

30.0 

30.0 

1 

175 

3 

0 

42 

50.2 

10.4 

.9 

32.0 

30.0 

1 

175 

3 

0 

42 

50.2 

11.4 

.9 

29.0 

29.0 

1 

175 

3 

0 

42 

50.2 

12.4 

.9 

24.0 

29.0 

1 

175 

3 

0 

42 

50.2 

13.4 

1.0 

23.0 

29.0 

1 

175 

3 

0 

42 

50.2 

14.4 

1.0 

26.0 

30.0 

1 

175 

3 

0 

42 

50.2 

15.4 

1.0 

25.0 

31.0 

1 

175 

3 

0 

42 

50.2 

16.4 

1.0 

27.0 

31.0 

1 

175 

3 

0 

42 

50.2 

17.4 

1.0 

36.0 

33.0 

1 

175 

3 

0 

42 

50.2 

18.4 

.9 

43.0 

35.0 

1 

175 

3 

0 

42 

50.2 

19.4 

.9 

44.0 

37.0 

1 

175 

4 

0 

40 

50.2 

7.4 

1.1 

17.0 

32.0 

1 

175 

4 

0 

40 

50.2 

8.4 

1.0 

16.0 

31.0 

1 

175 

4 

0 

40 

50.2 

9.4 

.9 

20.0 

30.0 

1 

175 

4 

0 

40 

50.2 

10.4 

.9 

24.0 

31.0 

1 

175 

4 

0 

40 

50.2 

11.4 

.9 

23.0 

30.0 

1 

175 

4 

0 

40 

50.2 

12.4 

.9 

23.0 

30.0 

1 

175 

4 

0 

40 

50.2 

13.4 

.9 

24.0 

31.0 

1 

175 

4 

0 

40 

50.2 

14.4 

.9 

24.0 

32.0 

1 

175 

4 

0 

40 

50.2 

15.4 

.9 

25.0 

33.0 

1 

175 

4 

0 

40 

50.2 

16.4 

.9 

30.0 

34.0 

1 

175 

4 

0 

40 

50.2 

17.4 

.9 

34.0 

36.0 

1 

175 

4 

0 

40 

50.2 

18.4 

.9 

37.0 

37.0 

1 

175 

4 

0 

40 

50.2 

19.4 

1.0 

35.0 

38.0 

1 

175 

4 

0 

38 

50.2 

6.4 

1.3 

27.0 

28.0 

1 

175 

4 

0 

38 

50.2 

7.4 

1.0 

23.0 

36.0 

1 

175 

4 

0 

38 

50.2 

8.4 

.9 

26.0 

35.0 

1 

175 

4 

0 

38 

50.2 

9.4 

.8 

33.0 

33.0 

1 

175 

4 

0 

38 

50.2 

10.4 

.8 

38.0 

34.0 

S8 


e.  Output  data  were  then  reformatted  and  moved  into  correctly 
named  files  and  directories  to  be  CCIR  compatible.  This  was 
accomplished  as  follows: 

1.  The  format  of  the  CCIR  compatible  field  strength  data 
base  file  for  use  in  the  statistical  filter  program  is  shown  in 
Table  3.  For  presentation  purposes  here,  the  format  in  the  table 
has  been  modified  to  include  column  descriptions  in  brackets. 
Table  3  is  a  condensed  form  of  the  complete  CCIR  data  base.  A 
sample  of  the  complete  CCIR  database  is  shown  in  Table  4  (format 
modified  with  bracketed  descriptors  for  clarity).  Path 
information  for  the  CCIR  database  is  given  in  Appendix  B. 

Table  3  also  includes,  in  the  last  column,  the  predicted 

value  of  field  strength  from  a  previous  version  of  the  FS  module. 
That  data  was  ignored  and  overwritten  with  the  new  predicted 
values,  as  described  below. 

2.  A  problem  with  the  format  of  Table  3  is  that  is  it  does 

not  list  a  unique  identifier,  such  as  universal  time  (UT).  A 
program  was  written  (Appendix  C)  to  merge  the  UT  and  other 

selected  parameters  with  the  data  in  Table  3  and  also  add  the 
predicted  FS  values  from  the  new  model  run.  A  sample  of  the  end 
result  is  shown  in  Table  5.  (Bracketed  descriptors  are  included 
for  clarity.) 

f.  Before  running  the  statistical  data  filter  (TURBO  program) 

which  generally  takes  40  to  50  hours  to  run,  various  manual 
checks  were  performed  to  reasonably  assure  that  the  run  would 

provide  useful  results.  These  checks  consisted'  of  spot 
comparisons  of  field  strength  outputs  from  the  modified  FS  module 
with  comparable  field  strength  data  in  the  CCIR  data  base.  If 

these  checks  showed  reasonable  comparisons  —  most  within  about 
10  to  20X  or  less  —  the  statistical  run  was  performed. 
Otherwise,  the  module  was  reviewed  to  determine  the  problem.  The 
problems  fell  into  various  categories,  including  logical  errors, 
incorrect  formats,  and  roundoff  errors.  Once  these  errors  were 
isolated  and  corrected,  the  statistical  run  was  performed. 

g.  The  statistical  run  consisted  of  running  the  statistical  data 
filter  program  using  the  formatted  data  as  described  above.  This 
run  took  approximately  40  to  50  hours  and  produced  line  printer 
listing  which  could  then  be  compared  visually.  To  aid  in  data 
analysis,  a  program  was  written  (mtable,  included  in  Appendix  D). 
This  program  reformatted  the  outputs  from  both  the  original 
hfbc84  runs  and  the  new  LTLFLD  outputs  into  summary  table 
formats.  An  analysis  of  the  these  tables  is  included  in  3.2.3. 
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Table  5.  Data  File  merg.hfbc  Sample 


Feb  07  21:28  1988  merg.hfbc 


[Data  base  code]  C101=HFBC84  (used  for  LTLFLD  model)] 
[Number  of  data  base  fields] 

[Total  data  base  records] 

[Data  base  title] 


12277  HFBC84  Field  Strenath  model 


[Path  ID] 

-  [Circuit  Length] 

[Season:  l=Winter,  2=Spring,  3=Summer,  4=Fall] 
ion:  0=other,  l=North/South,  2=East/Uest] 

•[Smoothed  SSN] 

-[Mid  path  latitude] 

-[Mid  path  local  time] 

[Operational  Freq/MUF  ratio] 
[—[Observed  Field  Strength] 

[Predicted  Field  Strength] 
[Universal  time] 

I —  [Month] 
l~[Year] 

[Operational  Freq/Operational 
[Monthly  Mean  SSN]. 


t 


MUFD 


1 

175 

3 

0 

42 

50.2 

G.4 

1.2 

23.0 

25.0 

6 

6.1 

8 

1984 

1.1 

25.0 

1 

175 

3 

0 

42 

50.2 

7.4 

1.0 

22.0 

27.0 

7 

6.1 

8 

1984 

1.1 

25.0 

1 

175 

3 

0 

42 

50.2 

8.4 

1.0 

22.0 

27.0 

8 

6.1 

8 

1984 

l.O 

25.0 

1 

175 

3 

0 

42 

50.2 

9.4 

.9 

30.0 

27.0 

9 

6. 1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

10.4 

.9 

32.0 

30.0 

10 

6.1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

11.4 

.9 

29.0 

30.0 

11 

6.1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

12.4 

.9 

24.0 

30.0 

12 

6.1 

8 

1984 

l.O 

25.0 

1 

175 

3 

0 

42 

50.2 

13.4 

1.0 

23.0 

27.0 

13 

6. 1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

14.4 

1.0 

26.0 

27.0 

14 

6.1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

15.4 

1.0 

25.0 

27.0 

15 

6.1 

8 

1984 

1.0 

25.0 

1 

175 

3 

0 

42 

50.2 

16.4 

1.0 

27.0 

27.0 

16 

6.1 

8 

1984 

1. 1 

25.0 

1 

175 

3 

0 

42 

50.2 

17.4 

1.0 

36.0 

26.0 

17 

6.1 

8 

1984 

1.1 

25.0 

1 

175 

3 

0 

42 

50.2 

18.4 

.9 

43.0 

24.0 

18 

6.1 

8 

1984 

1.2 

25.0 

1 

175 

3 

0 

42 

50.2 

19.4 

.9 

44.0 

19.0 

19 

6.1 

8 

1984 

1.3 

25.0 

1 

175 

4 

0 

40 

50.2 

7.4 

1.1 

17.0 

26.0 

7 

6.1 

9 

1984 

1.1 

15.0 

1 

175 

4 

0 

40 

50.2 

8.4 

1.0 

16.0 

27.0 

8 

6.1 

9 

1984 

1.1 

15.0 

1 

175 

4 

0 

40 

50.2 

9.4 

.9 

20.0 

27.0 

9 

6.1 

9 

1984 

1.0 

15.0 

1 

175 

4 

0 

40 

50.2 

10.4 

.9 

24.0 

27.0 

10 

6.1 

9 

1984 

1.0 

15.0 

1 

175 

4 

0 

40 

50.2 

11.4 

.9 

23.0 

26.0 

11 

6. 1 

9 

1984 

1.0 

15.0 

1 

175 

4 

0 

40 

50.2 

12.4 

.9 

23.0 

26.0 

12 

6.1 

9 

1984 

1.0 

15.0 

1 

175 

4 

0 

40 

50.2 

13.4 

.9 

24.0 

25.0 

13 

6.1 

9 

1984 

l.l 

15.0 

1 

175 

4 

0 

40 

50.2 

14.4 

.9 

24.0 

25.0 

14 

6.1 

9 

1984 

1.1 

15.0 
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3.2.3  Results  and  Conclusions 


The  results  of  the  PS  module  tests  are  shown  in  Appendix  E, 
tables  E-1  through  E-25.  Table  E~1  shows  the  overall  results  for 
all  81  paths  at  all  frequencies  recorded  in  the  database.  The 
results  are  satisfying  and  show  the  new  field  strength  method 
(LTFLD)  to  have  much  lower  average  residuals  than  the  old  method 
(DM6LT).  The  RMS  residual  is  also  better  by  at  least  a  factor  of 
two  to  one.  The  correlation  coefficient  is  better  but  the 
improvement  is  not  as  pronounced.  This  could  be  caused  by  a  time 
lag  in  the  prediction  relative  to  the  observed  data;  however  this 
has  not  been  verified.  The  most  pronounced  improvement  seems  to 
be  in  the  over-the-MUF  calculations  which  show  a  much  lower 
average  residual  of  approximately  one-fifth  and  an  RMS  residual 
of  less  than  one-half. 

These  results  look  encouraging;  however,  when  one  inspects 
Table  E-2  it  can  be  seen  that  the  shape  of  the  distribution  of 
losses  above  the  MUF  is  not  that  of  a  normal  curve,  but  possibly 
an  exponential  curve  as  suggested  by  Oamboldt  and  Sussman.  The 
scatter  fields  (well  above  1.4  times  the  MUF)  are  stronger  than 
predicted.  A  revised  method  for  over-the-MUF  computations  has 
been  postulated  but  time  and  resources  did  not  allow  checks  with 
the  data  base. 

Inspection  of  Tables  E-ll  through  E-13,  which  are  a  function 
of  path  length,  shows  that  an  improvement  has  been  made  for  the 
very  short  paths  (less  than  1000  km)  with  an  average  residual  of 
about  one-fifth  the  older  value.  This  was  4  major  concern  of  the 
project.  The  RMS  residuals  show  good  improvement;  however, the 
correlation  coefficient,  while  showing  good  improvement  in 
general,  looks  very  bad  for  paths  between  2000  and  3000  km.  It  is 
suspected  that  paths  of  this  length  should  be  easily  predictable. 
Further  inspection  shows  that  only  1.27.  of  the  data  base  was  used 
to  determine  these  values  of  comparison.  Also,  this  distance  is 
the  "crossover"  region  for  the  one-hop  E-layer  mode  and  the  one- 
hop  F2-layer  or  2-hop  E-layer  mode.  These  modes  change  from  day- 
to-day  due  to  changing  ionospheric  heights  and  critical 
frequencies.  This  may  help  to  explain  the  low  correlation 

coefficient,  but  this  has  not  been  analyzed  further. 

A  very  similar  trend  is  shown  between  4000-5000  km  ranges, 
which  is  the  "crossover"  point  for  one-hop  and  two-hop  f 2-layer 
propagation.  The  antenna  patterns  are  sensitive  to  these 

changing  modes  due  to  vertical  takeoff  angle  associated  with  the 
modes  available  day-to-day. 

Only  0.9)1  of  the  data  was  used  in  the  analyses  for  the  above 
range  of  4000-5000  km  suggesting  that  more  data  would  be  useful 
to  make  a  detailed  analysis  of  the  trend. 

The  over-the-MUF  calculation  shown  in  Table  E-12  as  a 
function  of  path  length  shows  large  average  residuals  for  the 
same  distances.  This  suqqests  the  same  conclusion  as  above 
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because  an  error  in  calculating  the  MUP  for  the  path  as  a  one  or 
two  hop  will  cause  a  more  than  usual  error  at  these  distances. 
The  data  tends  to  suggest  the  MUF  was  consistently  lower  than 
that  of  the  observed  data.  In  other  words,  over-the-MUF  losses 
were  calculated  when  in  fact  the  path  observed  was  operating 
below-the-MUF.  This  has  not  been  verified  and  it  is  doubtful 
that  data  is  available  to  verify  the  assumption. 

Another  possibility  is  the  data  recorded  may  be  only  on 
those  days  that  the  signal  was  readable  above  the  noise;  hence 
below  the  MUF  for  the  path.  If  all  30  days  at  a  given  hour  are 
not  present  in  the  day  it  is  difficult  to  tell  what  part  of  the 
distribution  it  represents.  There  is  no  intent  to  blame  the  data 
for  bad  predictions,  but  rather  to  cite  things  that  can  be  looked 
at  to  improve  the  overall  predicition  process. 

The  large  errors  occurring  when  very  few  hours  were  recorded 
may  be  due  to  an  insufficiently  large  sample  size  to  produce 
meaniful  statistics.  Although  included  here  for  completeness, 
they  are  bothersome.  It  is  suggested  that  these  data  be 

eliminated  for  analysis  purposes  that  do  not  a  least  come  close 
to  the  hourly  median  of  the  monthly  median  values,  even  though 
this  violates  the  desire  for  a  large  data  base.  The  smaller  data 
base  may  be  of  more  value  in  updating  the  prediction  methods. 

The  model  for  field  strength  predicitions  within  HFTDA  is 
hic^iily  modular  and  should  allow  easily  checking  different 
algorithms  for  all  the  losses  associated  with  the  system  loss  and 
field  strength  calculations. 

3.3  Recommendations 

The  availability  of  a  highly  modular  field  strength  model 
should  facilitate  checking  many  options,  specific  models  for 
given  parameters,  and  results  of  these  calculations.  We  suggest 
checking  any  changes  with  limited  circuits  from  the  CCIR 
database.  It  would  be  preferrable  to  use  those  circuits  with 
many  hours  of  data  even  if  some  path  lengths  are  avoided  in  the 
initial  checks  of  any  parameter. 

Worthwhile  items  to  be  revised  and  checked  with  observations 

are: 

1. )  Shultz-Gal let  absorption  equation 

2. )  Dambolt  over-the-MUF  exponential  distribution 

3. )  Revised  JTAC  model  for  scatter  fields  above  about  1.4«MUF 

4. )  Daytime  low-frequency  specular  reflection  algorithhm 

5. )  F2-layer  height  calculation 

6. )  High-latitude  description  of  losses 

7. )  Rise-time  of  the  MUF  of  MUF-85  near  mid-day 

8. )  E-layer  algorithm  for  MUF 
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These  are  not  necessarily  in  the  order  of  priority.  They 
represent  potential  areas  of  model  enhancement  which  could  lead 
to  expected  improvements  in  the  predictions. 

These  changes  should  be  made  to  see  which  algorithm  shows 
the  best  fit  with  the  overall  data,  to  make  HFTDA  a  true 
predictor  of  field-strengths.  Care  should  always  be  taken  to 
avoid  using  a  limited  amount  of  data  taken  in  a  specific  location 
to  draw  conclusions  aside  from  the  physics  of  the  matter.  A 
curve  fitting  process  which  is  not  founded  on  long-term  emperical 
data  and  physics  could  be  non-productive  while  showing  very  good 
statistical  data  of  which  many  transmitting  and  receiving 
variables  are  not  known;  especially  as  a  function  of  time. 

Since  the  MUF  and  vertical  arrival  angle  were  not  recorded 
at  the  receiver  in  the  CCIR  database,  the  power  and  elevation 
angle  gain  at  the  transmitting  antenna  are  of  questionable 
validity  over  a  period  of  time.  It  is  recommended  that  the  above 
remarks  be  kept  in  mind  when  deducing  the  validity  of  any  checks 
against  data  'to  determine  the  accuracy  of  HF  predictions  for 
their  varied  purposes. 
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APPENDIX  A.  FLOWCHART  SYMBOLOGY  AND  CONVENTIONS 


Figure  A-1  shows  the  flowchart  symbology.  In  general,  flowcharts 
use  FORTRAN  conventions  to  represent  discrete  programming  steps. 
Assignments  are  indicated  with  an  sign,  and  mean  the  right 

value  is  assigned  to  the  left  value.  The  right  value  is 
unchaged.  Thus, 

A  =  B 

means 


B  ->  A 

FORTRAN  arithmetic  symbols  are  used.  Thus, 

♦  -  multiplication 

+  -  addition 

-  subtraction 
♦♦  -  exponentiation 

/  -  division 

In  addition,  the  following  functions  or  their  equivalents 
are  assumed  to  be  available  or  easily  programmed,  and  so  are  not 
separately  charted. 

LOGIO  arithmetic  LOG  function,  base  10 

INT  floating  point  to  integer  type  conversion  function 

MAX, MIN  return  maximum  or  minimum  value  from  a  list 

COS,  SIN,  TAN,  SEC  the  trigonometric  functions 

ATAN 
SORT 
ABS 


the  inverse  tangent  trigonometric  functions 
the  square  root  function 
the  absolute  value  function 
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THE  FOLLOWING  SYMBOLS  HAVE  BEEN  USED  IN  THESE  FLOWCHARTS. 


CZZ) 


ENTRY  OR  EXIT  POINT 
USUALLY  LABELLED  WITH 
ROUTINE  NAME  OR  RETURN 


PROCESS 


CONDITIONAL  BRANCH 


I/O.  MASS  STORAGE 


CONNECTOR.  ON  OR  OFF  PAGE 


FOR  N  -  N1  TO  NN...NEXT 
STRUCTURE,  OR  DO  LOOP. 


COMMENTS 


Figure  A-1  Flowchart  Symbology 
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APPENDIX  B 

PATH  INFQRHATIQN  FOR  THE  CCIR  DATABASE 


Path  Information  for  the  HFBC84  Field  Strength  Database 
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Path  information  for  the  HFBC84  Field  Strength  Database  (cxint) 
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APPENDIX  C 
DATA  MERGE  PROGRAM 


Feb  07  21:09  1988  rorgfs.f 


program  mrgfs 
c 

c  This  program  merges  Field  Strength  data  generated  by  the  TDA  program 
c  FSTREN  into  the  field  strength  CCIR  database  'fsdb.hfbc’,  based  on  the 
c  universal  time,  frequency  and  SSN,  as  indexed  by  control  file  'fsccir4’ 
c  and  saves  the  final  merged  output  into  'merg.hfbc’ 
c 

implicit  none 
integer  i dummy 
c 

c  fsccirA  declarations  (note:  ut:  range=l-24) 

real  ut , month , year , mmssn , pat hi n , r 1 at , r 1  on , mi dpt , 1 1  at , 1 1 on , 

*  dummy, opfreq, ofqluf , lufmuf , fs,ofqomf ,rmssn, luf ,ofqmuf ,opmuf 
integer  pathid,rec4 

c 

c  fsdb.hfbc  header  declarations 

integer  is,nflds,nrecs,recdb 
characterise  dbname 
character *80  fmt 
c 

c  fsdb.hfbc  data  record  declarations 

real  mlat, It, fmuf,obs,predl,sssn 
integer  id, len,sea,ort,ssn 
c 

c  tstnew  (predicted  model)  declarations  (note:  Imt:  range=0-23) 

integer  ssnp,monp,yearp, loop, itp,recp 
real  mmssnp,lmt, fst, frq,luf2,ffluf2,omuf2,utp, fstO 

c 

c  merg.hfbc  declarations 
integer  reem 
character*80  fmtnew 
data  fmtnew 

*  /’ (i2,i6,2i2,i4,f7.1,f5.1,f4.1,2f6.1  , i3, f5. 1, i3, i5,2f6. 1) ' / 
c 

c  initialize  counters 
rec4=0 
recdb=0 
recp=0 
c 

c  open  all  files  used  by  this  program 

opendO,  file='  fsdb.hfbc'  ,status='old'  ,err=9100) 
opendl,  file='tstnew'  ,status='old'  ,err=9110) 
open(12, file®' fsccir4'  ,status='old' ,err=9120) 
open(13, file='merg.hfbc' , status='new' ,err=9130) 
write(*,986) 

986  format ('  files  opened') 
c 
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c  copy  over  database  header 

read(10,80,err=9100)  is,nflds,nrecs,dbnamef  fmt 
80  format (215, 18, 2x, a/a) 

wrlte(*,984)  _ 

984  formate*  header  read') 

wrlte(13,80)  ld,nflds,nrecs,dbname, fmtnew 
recm=l 

wrlte(*,988) 

988  formate  header  written’) 
c 

c  read  In  Initial  observed  database  record 

r ead ( 1 2 , 1 1 0 , er r  =9 1 20 , end=9 1 20 )  ut , month , year , mmssn , pat h 1 n , 

*  rlat,rlon,mldpt,tlat ,tlon,dummy,opfreq,ofqluf , lufmuf , fs, 

*  ofqomf , rmssn, luf ,ofqmuf ,opmuf ,pathld 
rec4=3 

no  format(2(7fl0.0,/),6fl0.0,7x,13) 

read (10, fmt ,err=9100,end=9100)  Id, len, sea, ort , ssn,mlat , It, 

*  fmuf , obs,predl 
recdb=2 
sssn=ssn 
wrlte(*,989) 

989  formate  Initial  db  record  read  In’) 
c 

c  make  sure  field  strength  databases  are  In  sync... 

If  (rmssn.ne.sssn.or. fs.ne.obs.or.pathld.ne.ld)  then 
wrlte(*,140)  rec4,recdb 

140  formate  mismatch  at  ccir4  db  (’,14,’)  and  hfbc  (’, 

*  14,’)’) 
endl  f 

c 

c  read  In  values  from  predicted  model  run 
145  readdl,  150,err=9110,end=5000)  ssnp,monp,yearp, dummy 

ISO  format (39x, 13,5x, 12, lx, 14/i2x, fS.O) 

recp=r ecp+1 
mmssnp=ssnp 
c 

c  loop  through  all  24  entries  searching  for  a  Universal  Time  match 
c  (note  that  for  loop=0,  utp=0,  but  fscclr4  doesn’t  have  ut=0,  but  Instead 
c  has  ut=24. . .therefore,  on  loop  #  ’O’,  save  fs,  and  for  loop  #  ’24’  use 

c  that  value  as  If  utp=24!  If  you  thought  that  was  complicated, 
c  you  should  see  how  long  It  took  to  figure  this  out!!) 
do  300  loop=0,24 
If  (loop.eq.24)  then 
utp=24.0 
fst=fst0 

else 

readdl,  160,err=9110,end=9110)  itp,lmt,  fst,  frq, 

*  luf2,ffluf2,omuf2 

160  format (1 10, £f 10.5) 

recp*recp+l 

utp^ltp 

If  (Itp.eq.O)  fst0=fst 
endl  f 
c 

c-a. 


c  check  for  a  match  to  the  fs  database 

if  (utp.eq.ut.and. frq.eq.opfreq.and.mmssnp.eq.mmssn)  then 
c 

c  found  a  match  -  create  updated  entry  in  merged  database 

wr ite( 13, fmtnew) id, len,sea,ort ,ssn,mlat , 1 t , fmuf ,obs, 

*  fst,utp, op freq, month, year , ofqmuf ,mmssn 
recm=recm+l 

c 

c  read  in  next  database  entry 

r  ead  < 1 2 , 1 1 0 , er  r  =9 1 20 , end=5000 )  ut , mont  h , year , mmssn , 

*  pathln,rlat ,rlon,mi dpt, tlat,t Ion, dummy, op freq, ofqluf, 

*  lufmuf, fs, ofqomf,rmssn,luf, ofqmuf, opmuffpathid 
rec4=rec4+3 

r  ead ( 1 0 , f mt , er  r =9 1 00 , end=5000 )  i d , 1 en , sea , or  t , ssn , ml  at , 

*  It, fmuf ,obs,predl 
recdb=recdb+l 
sssn=ssn 

idummy=mod (recdb, 100) 

if  (idummy.eq.O)  write(*,987)  recdb, rec4 
987  formate'  recdb=’,i6,'  (12279),  rec4=',i6,'  (36831),  ') 

c 

c  are  field  strength  databases  still  in  sync?.... 

if  (rmssn.ne.sssn.or. fs.ne.obs.or .pathid.ne. id)  then 
write(*,140)  rec4, recdb 
endi  f 

endi  f 

300  continue 

c 

c  all  entries  from  current  configuration  of  predicted  model  are 
c  complete  -  go  read  in  another  predicted  model  data  set 
go  to  145 
c 

c  eof  encountered  on  any  cf  the  input  files  has  occured 

5000  write(*,5010)  recdb, recp,rec4,recm 

5010  format (lx, i5, ’  recs  read  from  fsdb.hfbc'/ 

*  lx,i5,’  recs  read  from  tstnew’/ 

*  lx,i5,'  recs  read  from  fsccir4'/ 

*  lx,i5,’  recs  written  to  merg.hfbc') 

close( 10) 

c 1 ose (11) 
close( 12) 
close( 13) 
go  to  9900 
c 

c  disk  errors 


9100 

write(*,9105) 

9105 

format (’  disk  error 
go  to  5000 

on  fsdb.hfbc’/) 

c 

9110 

write(*,9115) 

9115 

format ('  disk  error 
go  to  5000 

on  tstnew'/) 

c 
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9120  write(*,9125) 

9125  formate’  disk  error  on  fsccir4’/) 

go  to  5000 
c 

9130  write(*,9135) 

9135  formate'  disk  error  on  merg.hfbc’/) 

go  to  5000 
c 

9900  end 


APPENDIX  D 
Statistical  Table 
Generation  Program 


program  mtabl 
c 

c  This  program  reads  thru  a  "DSOxx"  statistical  output  file  generated 
by 

c  the  TURBO  program  and  creates  an  ascii  table  for  use-±ft- subsequent 
c  report  generations.  Two  input  lines  are  required: 

c  "DSOxx"  filename  (80  char)  and  output  'table'  filename  (80  chars) 

c  The  output  file  (table)  is  in  the  same  format  as  that  generated  in 
the 

c  original  1984  report.  The  LTLFLD  columns  are  filled  in  by  reads  to 
the 

c  "DSOxx"  file  located  in  you  current  directory  and  the  DAMBLT  columns 
are 

c  filled  in  by  reads  to  the  "DSOxx"  file  located  in  the  original  field 
c  strength  directory  ' /bigdisk/datscr/fldstr/results/ "DSOxx" ' . 
c 

implicit  none 
real  totpop 

parameter  (totpop=12277 . 0 ) 

integer  ncmax 
parameter  (ncmax=8) 

integer  mxline 
parameter  (mxline=30) 

integer  ncvau:  (ncmax)  ,popsz,idtBTimy  (10)  ,i,nconds,  line  ' 
integer  zcvar (ncmax) , zopsz 
character*!  rangesep 

character*2  codes (ncmax) , zodes (ncmax) ,  ch 

real  cval (ncmax ) , AvgRes , RmsRes , dummy ,  AvgRelRs ,  RmsRelRs ,  AvAbRlRs 
real  zval (ncmax) ,  zvgres, zmsres,  zvgrelrs, zmsrelrs, zvabrlrs 

real  dummy3 (3) , Corr, PcPop, zorr, zcpop 
character*65  infyl, oufyl,apfyl 
character* 31  apbeg 

data  apbeg/ '  /bigdisk/datscr/fldstr/resxjlts/ '  / 
write (6, 140) 

140  formate  Enter  input  filename  <rtn>  output  (table)  filename') 
read(5, 160,err=900)  infyl, oufyl 
160  format (a/a) 

apfyl=apbeg  //  infyl 
write (6, 200)  infyl, oufyl, apfyl 
200  formate  Stat  file  (input)  is  :',A  / 

*  '  Table  file  (output)  is  :',A  / 

*  '  APES  (NOSC)  file  is  :',A) 
c 

c  open  all  used  files 

open (10, file=infyl, status='old' ,err=9100) 
open (11, file=oufyl, status='new' ,err=9110) 
open (12, file=apfyl, status='old' ,err=9120) 
line=mxline 
c 

c  read  in  a  logical  data  'record'  (5  physical  records)  from  "DSOxx" 
file 

300  readdO,  340,err*9100,end=800)  (ncvar  (i) ,  codes  (i) ,  cval  (i) ,  i=l,  8) 
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340  fonnat(8(i2,a2,fl0.2,lx)) 


read (10,  380,  e  r=9100)  AvgRes,RmsRes,  dummy,  AvgRelRs,RinsRelRs, 

*  dummy,  AvAbRlRs, dummy3,Corr,popsz,  dxrnmy3, dummyS, 

*  dummy3, idummy 

380  format (6gl0 . 4e2/5gl0. 4e2, il0/6gl0. 4e2/2gl0 . 4e2, f 6. 3, 16, 14,812) 

c 

c  count  number  of  conditions  satisfied  to  determine  printout  format 
nconds=0 
do  l=ncmax,l,-l 
If  (codes (1) . ne. '  ’)  then 

nconds=l 
go  to  460 
endlf 

enddo 

c 

c  read  In  a  logical  data  'record*  from  original  APES  NOSC  file 
460  read(12, 340, err=9120)  (zcvar (1) , zodes (1) , zval (i) , 1=1, 8) 

read (12, 380, err=9120)  zvgres, zmsres, dummy, zvgrelrs, zmsrelrs, 

*  dummy,  zvabrlrs,  duinmy3,  zorr,  zopsz,  dummy3,  dijmmy3, 

*  dummy  3 ,  Idimmy 
c 

c  calculate  percent  of  total  population 
dummy=popsz 

PcPop= (dummy/totpop) *100 . 0 
dummy=zopsz 

-cpop= (dummy/totpop) *100 . 0 
llne=llne+l 

If  (llne.le.mxllne)  go  to  540 
c 

c  write  out  table  header  first 
write (11, 500,err=9110) 

500  format  (lhl/2x, 'Conditions' ,  4x, 'Avg  Residual ',  5x, 'Rms  Resldioal ' , 

*  5x, 'Avg  Rel  Res  ',5x, 'Rms  Rel  Res  ',3x, 'A^'^g  Abs  Rel  Res  ', 

*  2x, 'Correlation  ',3x, '%  of  Total' 

*  /  llx,  5(4x, 'LtFld',3x, 'DmBlt'),  2 (4x, 'LFld' , 3x, 'DBlt ' ) 

*  /  128('»')  ) 
llne=l 

c 

c  determine  table  printout  format 
c 

c  None  Single  Multiple 

540  go  to  (  580,  660,  740  )  nconds+1 

If  (nconds.gt.O)  go  to  740 
c 

c  no  conditions  found  (assume  OverAll  run) 

580  write (11, 620, err=9110)  AvgRes, zvgres, RmsRes, zmsres, AvgRelRs, 

*  zvgrelrs, RmsRelRs, zmsrelrs, AvAbRlRs, zvabrlrs, Corr, zorr, 

*  PcPop, zcpop 

620  format (12x,5('  I  ',f5.1,'  /  ',f5.1,lx), 

*  '  I  ',f4.2,'  /  ',f4.2,lx,'  I  •,f4.1,'  /  ',f4.1) 
go  to  300 

c 

c  one  condition  found  -  ■'ust  print  it  out 
660  ch='=  ' 
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if  (codes  (nconds)  .eq. 'LT' )  ch=''<  ' 
if  (codes (nconds) .eq. 'LE' )  ch='<=' 
if  (codes (nconds) .eq. 'GT' )  ch®'>  ' 
if  (codes  (nconds)  .eq. 'GE' )  ch*='>=' 
if  (codes (nconds) .eq. 'NE* )  ch»’<>' 

write  (11,  700,err=9110)  ch,cval  (nconds)  ,AvgRes,  zvgres,RinsRes, 

*  zmsres,AvgRelRs, zvgrelrs,RmsRelRs, zmsrelrS/AvAbRlRs, 

*  zvabrlrs, Corr, zorr, PcPop, zcpop 

700  format (lx,a2,lx,f 8.2,5 (•  1  ',f5.1, *  /  ',f5.1,lx), 

*  '  I  ',f4.2,'  /  ',f4.2,lx,'  I  •,f4.1,'  /  ',f4.1) 
go  to  300 

c 

c  multiple  conditions  found  -  -  assume  last  two  are  to  be  printed 
740  if  (codes (nconds) .eq. 'EQ' )  go  to  660 
rangesep='-' 

if  (cval (nconds-1) .ge. 1000. 0. or. cval (nconds) .ge. 1000. 0)  then 
cval (nconds-1) =cval (nconds-1) /lOOO.O 
cval (nconds) =cval (nconds) /I 000. 0 
rangesep='=' 
endif 

write (11, 780,err=9110)  cval (nconds-1) ,rangesep, cval (nconds) , 

*  AvgRes, zvgres,  RmsRes, zmsres,  AvgRelRs, zvgrelrs, RmsRelRs, 

*  zmsrelrs, AvAbRlRs, zvabrlrs, Corr, zorr, PcPop, zcpop 
780  format (lx,f 5.1, al,f 5.1,5 (’  |  ',f5.1,'  /  ',f5.1,lx), 

'  I  ',f4.2,'  /  ',f4.2,lx,'  I  •,f4.1,'  /  ',f4.1) 
go  to  300 
c 

c  EOF  on  data 
800  write (11, 820) 

820  format (Ihl) 
go  to  9900 
c 

c  error  on  command  input 
900  write  (6,  910) 

910  formate  Error  reading  commands  (filenames)') 
go  to  9900 
c 

c  error  on  input  (DSOxx)  file 
9100  write (6, 9105)  infyl 

9105  formate  Error  reading  data  from  DSOxx  file  ',A) 
go  to  9900 
c 

c  error  on  output  (table)  file 
9110  write (6, 9115)  oufyl 

9115  formate  Error  writing  to  table  file  ',A) 
go  to  9900 
c 

c  error  on  orig  APES  NOSC  file 
9120  write (6, 9125)  apfyl 

9125  formate  Error  reading  data  from  orig.  NOSC  "DSOxx"  file  ',a) 
c 

c  close  all  files  and  exit 
9900  close (10) 
close (11) 
end 
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Field  Strength  Test  Results 
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Table  e-2.  Summary  of  results  as  a  function  of  Frequency  to  MUF  ratio  for  the  entire  data  base 


Conditions  Avg  Residual  Rms  Residual  Avg  Rel  Res  Rms  Rel  Res  Avg  Abs  Rel  Res  Correlation  %  of  Total 

LtFld  DmBlt  LtFld  DmBlt  LtFld  DmBlt  LtFld  DmBlt  LtFld  Dmblt  LFld  DBlt  LFld  DBlt 


mr^o>rof^»HfHCsjiHCN  v 


^rHoooi/)^in^r-a^fHiH^csjo^rooinrooorooOf-»i-»OiH^iHOiHfHr^a»ocMincvj 


ro  r**  ro 


rg  fH  eg 


^^>^CT^rgcDa^oc^c^ooo^^cD^cgegvDovDo^HO^lnoocgr“CO•Hcg^»HO<^^'O^^Of-l^oo 

^-HCDvfimr^l/)vvo^or«•r•^nvomooaovcooo^^coo^rgvom^c^^*lHr^mor*•vo^^r>oovoo 


^-^rgm^^^^01-Ha^c^f-^cglnoooo^Hcg«-lOf^lnoo^ocg«  o^r^oor^ooi-icDCDOvovooocgoo 
^ao<^fnco^nof^fnaoo^a^^fir-^ooooooo^oo^^-^^^gt  eg^ootnor-4r*>momro^fH4-^vD 


tn^m^-»^POvorolr>l-•»9'voc»mmoo<nc^Jvcovo^ovD^AfHv^cOfHf-HfHc^fnrgvcNO^rgu^^ 


eg  in  m  VP  rs  ^ 


m  «-g  lO 


r*  rH  ro 


ro^fn*-i»-iro^ooo^^<NVoegfnin^^roen^i/>eg4nvPi-<if>egr»^egr*i-gf-^egvoif>«Hfomvp 


CM  m  ^  min 


m  CD 


CM  m 


tnr^m«-HOOVDmr^avinr-invpm«-H^CMmegvm^vpi/>#H\PorgiOegrMroeMvpi/>vmr*^ 


^  m  o 

VP  ^ 


movpcMvnvpmvpcM^mm^inr»*fHvp^ovvommv^egmegmvpcot^fHmoovpvpvir>vp 

V  ^r*fHmcM  om^  o  O  m  oovm 

^eM»HCN<N*Hf-g  CMCM  ^^eg 


^^Ot-«r-ovpininaveMv©vegvncoegegi-40^0f*g  vo  in  ovpooocM*Hegmovp^^*-imv 


VO  i-g  <s  I  I 


VP  I  m 
I  I  I 


«-t^rgi-ir-f^^vnvpco^cgeMO^in^^cgeg^u^iHmvPt-iin^povp^iH»^rgvp^<nmavvp 


r*Ofnegmvoop^mr^r^«Hcginooooino^m^egtHO%mv'o^veg4negoin*Havavinvpi-g 

lnocDmooromr“CS|^oomo^oovpQo*-^lnegmcnvocMr*^omav^commQC>r•mvpvpvPGoeg 
•-ifH  fOfHCNromm  f-t  m  tH  eHegrg«H  csif-i  eg  »h  ih»hih  »h 


vvp(M,-Hpv^vpoegmvpp»avocMvpvpovp^o»HvpooDegiHr**omocDODi-4o«H'^ininfH 

roooolnmm^^o^moD^GO'0^vpeg^p^commr•o^eMvr^^OP*•^m^MOO^nvpop-o 


vpc^o^Olnp^a^o^vpcooooopo^HoovplnoOfHr•spoor•Olnvo^cNOcooo^^^vpcgvp^m 


egmooegOfHO%*Heg«-«egm^r-in^oegegpoegr*- 
eg  egeg  ih  <Hitegi  egegi-i  rM*H 

I  I 


m  ro  m  r*  eg  o 

1  fH  I  I  tH 


^  ro 
I  I  ^ 


eMODegegr•oo^m^noop-►m4nvp,-^voo^^eg•HmooooDa^mlnlnm^megoOlnvpr-ln 


inooeg^vp^mvpm^vpm 

II  fH  1  I 


m  I  o  ^  eg  vn  eg  fH  « 
fH  rH  I  eg  I 


'  CD  o  m  VP  eg  00 
fH  eg  iH 


comegvpvpiHvoo 


0000000000000300000000000000000000000000 

oooooooooooooooooooooooooooooooooooooooo 

fHCMf0^vnvpc^oochOfHegm^in»pr-oooo^egmvinvpp»*ooovo^egm^4nvpp«>oDOveg 

,-4rHfHfHfHfHrH^^^ogegegegrgegegegegegmmmr)fnmmmmro^ 


I  I  I  I  I  I  I  I  I  I 


I  I  I  I  I 

E-4 


I  I  I  I  I  I 


I  I  I  I  I  I  I  I 


Table  e-3  .  Summary  of  results  as  a  function  of  Circuit  ID  for  the  entire  data  base 
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Table  e-4  .  Summary  of  results  as  a  function  of  Circuit  ID  for  the  entire  data  base 
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Table  E-5.  Summary  of  results  as  a  function  of  Circuit  ID  for  Frequency-to-MUF  ratio  >  1.0 
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Table  e-6.  Summary  of  results  as  a  function  of  Circuit  ID  for  Frequency-to-MUF  ratio  <=  1.0 
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Summary  of  results  as  a  function  of  Circuit  ID  for  Frequency-to-MUF  ratio  <=  1.0 


Conditions  Avg  Residual  Rms  Residual  Avg  Rel  Res  Rms  Rel  Res  Avg  Abs  Rel  Res  Correlation  %  of  Total 
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Table  e-8.  Summary  of  results  as  a  function  of  Midpath  Latitude  for  entire  data  base 


Table  e-9.  Summary  of  results  as  a  function  of  Midpath  Latitude  for  Frequency-to-MUF  ratio  >  1.0 
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Table  e-io.  Summary  of  results  as  a  funciton  of  Midpath  Latitude  for  Frequency-to-MUF  ratio  <=  1.0 
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Table  E-ll.  Summary  of  results  as  a  function  of  Circuit  Length  for  entire  data  base 


E-13 


Table  E-12.  Summary  of  results  as  a  function  of  Circuit  Length  for  Frequency-to-MUF  ratio  >  1.0 
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TableE-13.  Summary  of  results  as  a  function  of  Circuit  Length  for  Frequency-to-MUF  ratio  <=  1.0 
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Table  E-14.  Summary  of  results  as  a  function  of  Local  Time  for  entire  data  base 
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Table  E-15.  Summary  of  results  as  a  function  of  Local  Time  for  Frequency-to-MUF  ratio  >  1.0 
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TablcE-ia  Summary  of  results  as  a  function  of  Local  Time  for  Frequency-to-MUF  ratio  <=  1.0 
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Table  E-17.  Summary  of  results  as  a  function  of  Orientation  for  entire  data  base 
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TableE-18.  Summary  of  results  as  a  function  of  Orientation  for  Frequency-to-MUF  ratio  >  1.0 
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Table  E-19.  Summary  of  results  as  a  function  of  Orientation  for  Frequency-to-MUF  ratio  <=  1.0 
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Table E-20.  Summary  of  results  as  a  function  of  Season  for  entire  data  base 
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Table E-2L  Summary  of  results  as  a  function  of  Season  for  Frequency-to-MUF  ratio  >1.0 
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TableE-2Z  Summary  of  results  as  a  function  of  Season  for  Frequence-to-MUF  ratio  <=  1.0 
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TableE-23.  Summary  of  results  as  a  function  of  Smoothed  Sunspot  Number  for  entire  data  base 
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Table  E-24.  Summary  of  results  as  a  function  of  Smoothed  Sunspot  Number  for  Frequency-to-MUF  ratio  >  1.0 
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TableE-25.  Summary  of  results  as  a  function  of  Smoothed  Sunspot  Number  for  Frequency-to-MUF  ratio  <=  1.0 
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